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ABSTRACT 

In this paper we evaluate the spectrum of the pregalactic density field on scales 
l/i -1 Mpc < r < 100/i _1 Mpc from a variety of astronomical data. We start with the 
APM data on the projected angular correlation function, w(9), in six narrow magnitude 
bins and check whether possible evolutionary effects can affect inversion of the w(6) 
data in terms of the underlying power spectrum. This is done by normalizing to the 
angular correlation function on small scales where the underlying 3-dimensional galaxy 
correlation function, £(r), is known. Using the APM data in narrow magnitude bins 
allows us to test the various fits to the APM data power spectrum more accurately. We 
find that for linear scales r > 10/i _1 Mpc the Baugh and Efstathiou (1993) spectrum of 
galaxy distribution gives the best fit to the data at all depths. Fitting power spectra of 
CDM models to the data at all depths requires f2/i = 0.2 if the primordial index n = 1 
and Qh = 0.3 if the spectrum is tilted with n = 0.7. Next we compare the peculiar 
velocity field predicted by the APM spectrum of galaxy (light) distribution with the 
actual velocity data. The two fields are consistent and the comparison suggests that 
the bias factor is scale independent with O a6 /6 ~(0.2-0.4). These steps enable us to 
fix the pregalactic mass density field on scales between 10 and ~ 100/i _1 Mpc. The 
next dataset we use to determine the pregalactic density field comes from the cluster 
correlation data. We calculate in detail the amplification of the cluster correlation 
function due to gravitational clustering and use the data on both the slope of the 
cluster correlation function and its amplitude-richness dependence. Cluster masses 
are normalized using the Coma cluster. We find that no CDM model can fit all the 
three datasets: APM data on w(9), the data on cluster correlation function, and the 
data on the latter's amplitude-richness dependence. Next we show that the data on 
the amplitude-richness dependence can be used directly to obtain the spectrum of the 
pregalactic density field. Applying the method to the data, we recover the density field 
on scales between 5 and 25/i _1 Mpc whose slope is in good agreement with the APM 
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data on the same scales. Requiring the two amplitudes to coincide, fixes the value of 
Q to be 0.3 in agreement with observations of the dynamics of the Coma cluster. We 
then use the data on high- z objects to constrain the small-scale part, (l-5)/i _1 Mpc of 
the pregalactic density field. We argue that the data at high redshifts require more 
power than given by CDM models normalized to the APM and cluster data. Then we 
reconstruct the pregalactic density field out of which modern-day galaxies have formed. 
We use the data on blue absolute luminosities, the fundamental plane relations and the 
latest X-ray data on the halo velocity dispersion. From this we recover the pregalactic 
density field on comoving scales between 1 and 5/i~ 1 Mpc which is in reasonable 
agreement with the simple power-law extrapolation from the larger scales. 

Subject headings: cosmology: theory — cosmology: dark matter — galaxies: formation 
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1. Introduction 

The origin of structure in the Universe is one of the most outstanding problems in modern 
cosmology. In the gravitational instability picture it is assumed that structures in the Universe 
formed by gravitational growth of small density fluctuations seeded at some early epoch of the 
Universe evolution. The COBE discovery of the microwave background anisotropies (Smoot et al 
1992, Bennett et al 1996) proved convincingly that density fluctuations were already present at 
z ~ 1000. It is then reasonable to assume that these were indeed the seeds of the density field 
that were to lead to the present day structures. Following the COBE discovery of the large-scale 
structure at z ~ 1000 it is therefore even more imperative to try to reconstruct the pregalactic 
density field on scales and at epochs inaccessible to the microwave background measurements. 

The spectrum of the density field determines the epoch and the order of galaxy and structure 
formation. Its normalization point is fixed by observations that show that fluctuations in the 
galaxy counts today at z = have unity amplitude on scale rg = 8/i _1 Mpc. Given this amplitude 
at r§, or mass scale ~ 1O 15 M , and the slope of the spectrum, one can predict when and on what 
scales the first objects in the Universe began to collapse. The converse is also true. Similarly, on 
scales > r§ where the density field is still in linear regime, the various properties of galaxy and 
cluster distribution can be used to determine the spectrum of the density field on scales currently 
inaccessible to space-borne microwave background measurements. 

On theoretical level, the density field is assumed to have been seeded at some early epoch 
in the evolution of the Universe. The most popular of such mechanisms involves inflationary 
cosmology with cold-dark-matter contributing most of the mass in the Universe. Various 
discussions have also concentrated on the density fluctuations seeded by topological defects, 
such as strings and textures, that can be produced naturally during early phase transitions. 
Another possibility is the empirical approach, not necessarily motivated by the current ideas from 
high-energy physics, whereby one deduces the properties of the early Universe from the present 
data, instead of "predicting" what the present-day Universe should be. 
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The density field at some early time can be characterised by both its spatial and Fourier 
components which are related via 5{x) = (2ir)~ 3 j 5^ exp(— ik ■ x)d 3 k. If one assumes the density 
field to be a random variable, one can describe its statistical properties via the moments of its 
probability distribution. Assuming an overall spherical symmetry the correlation function of the 
field is defined as the generalized scalar product £(r) = (5(x + r)5(x)}. The power spectrum is 
defined as P(k) = (|<5^| 2 ), where the average is performed over all phases. The correlation function 
and the power spectrum represent a pair of three dimensional Fourier transforms. In addition 
the mean square fluctuation over a sphere of volume V containing mass M is related to £(r) via: 
A 2 (M) = (( J- ' s ^ dV ^j j — I t(r) dV _ jf phases of 5^ are random, the distribution of the density 
field is Gaussian, and the correlation function (or its Fourier transform) uniquely describes all 
properties of the density distribution. The central theorem ensures that the Gaussian distribution 
is reached in most inflationary mechanisms for generating energy-density fluctuations. In models 
with topological defects, however, the density field is not Gaussian and, in addition to £(r), it is 
determined by higher moments of the probability distribution. 

Enough data have by now been accumulated over a large range of scales to strongly constrain 
the models, on the one hand, and on the other to attempt to determine the spectrum of the 
pregalactic density field independently of theoretical prejudices. This is the aim of this article 
where we analyse the various data to show that the data lead to a consistent and unique density 
field in the pregalactic Universe. 

When using the data derived from galaxy catalogs one determines the distribution of light, 
whereas in order to understand the physics of galaxy formation and early Universe one needs 
to determine the distribution of mass. Therefore, one has to introduce the biasing factor which 
determines how the two are related to each other. The most common assumption is that of linear 
basing, <5ii g ht oc <5 mass , although more complex and complicated biasing schemes can be considered. 
Therefore, before we proceed to the main results of the paper we address the reliability of the 
general assumption of light-tracing-mass by comparing the velocity field predicted by the galaxy 
correlation function data with peculiar velocities in the Great Attractor region. 
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The plan of the paper is as follows: In sec. 2 we briefly define the cold-dark-matter models 
which we will test against the data in the paper and discuss the microwave background (COBE 
DMR) constraints on the very large scale part of the density field where it is likely that the 
power spectrum has preserved its original slope. In Sec. 3 we analyse the data from the APM 
survey divided into narrow magnitude bins (Am ~ 0.5) which contain galaxies at the same 
evolutionary stages. In this way we establish whether significant galaxy evolution effects could 
affect determination of the spectrum of the density field from the projected APM data. Following 
this, we discuss the various fits to the APM data and conclude that they all give essentially the 
same field of light (galaxy) distribution on linear scales, r > r§. The data in the narrow magnitude 
bins, however, set somewhat stronger limits on CDM models than used before. In Sec. 4 we 
elaborate on a method to interrelate the data on the galaxy correlation function from both APM 
and CfA catalogs to that on peculiar velocities. This establishes if the light traces light and the 
constraints both datasets suggest for f2 and the bias factor. Thus we determine the spectrum of 
the pregalactic density field and its correlation function on scales (10-100)/t _1 Mpc. In Sec. 5 we 
calculate the amplification in the cluster correlation function due to gravitational clustering of 
density fluctuations. We then use the data on both the spatial slope of the cluster correlation 
function and on its richness-amplitude dependence in order to set further constraints on the CDM 
models. We devise a method to invert the richness - correlation amplitude dependence of clusters 
of galaxies in order to obtain directly the spectrum of the pregalactic density field on scales 
(5-20)fo _1 Mpc. In Sec. 6 we provide comparison between the density field deduced independently 
from the APM catalog and from the cluster correlation data in order to constrain by requiring 
that the two match on the same scales. In sec. 7 we discuss the limits on the small-scale part of the 
pregalactic density field from the data on the existence and ages of high-z galaxies and clusters. 
We also reconstruct the pregalactic density field from the data on the fundamental plane of 
modern-day (elliptical) galaxies, their absolute photometry and the dynamics of their haloes from 
the recent X-ray observations. At the end of this, we reconstruct in Sec. 8 the pregalactic density 
field over two decades in linear scale, ~(l-100)/i _1 Mpc. Conclusions are summarized in Section 9. 



- 7- 



2. The early Universe anzatz: microwave background and very large scales 

Because of its elegance and simplicity inflation is probably the most popular of the current 
theories for the origin of structure in the Universe. The density field in inflationary picture 
originates from quantum energy fluctuations generated during the slow roll-over of the inflaton 
field during this era. Exponential expansion of the underlying space-time during inflation ensures 
that the fluctuations are seeded by normal causal processes on scales well in excess of the current 
particle horizon of ~ 6000/i -1 Mpc. In its simplest and most natural form inflation leads to 
gaussian adiabatic density fluctuations with the initial power spectrum which is scale-free and of 
the Harrison-Zeldovich slope, Pi(k) oc k n with n = 1. In conjunction with the COBE observations 
inflationary models generally require flat Universe (Kashlinsky, Tkachev and Frieman 1994) 
while Big Bang nucleosynthesis implies a small baryon density (Pagel 1997, Walker et al 1991). 
Thus in order to reconcile inflationary prejudices with observations one has to postulate the 
existence of cold-dark-matter that was not directly coupled to the baryon-photon plasma in the 
prerecombination Universe (e.g. Blumenthal et al 1984 and references cited therein). 

Dynamical evolution of cold-dark-matter (CDM) in the early Universe is well understood 
by now (Peebles 1982, Vittorio and Silk 1985, Bond and Efstathiou 1985). In the CDM model 
one assumes that evolution during some early (inflationary) epoch resulted in adiabatic gaussian 
density fluctuations with initially scale-invariant power spectrum. Later the evolution of density 
fluctuations leads to modification of the power spectrum because of the different growth rates of 
sub- and super-horizon harmonics during the radiation-dominated era. Thus in CDM models the 
power spectrum of the density field at the epoch of recombination is given by: 

P{k) oc k n T 2 {k) (1) 

The transfer function, T(k), depends mainly on the size of the horizon scale at the matter-radiation 
equality, oc (Qh 2 )^ 1 , and is 

= ln(l + a fc) + aifc + 2 + 3 + (a4jfe)4] -i/4 (2) 
aok 

where a = 0.6ai = 0.14a 2 = 0.43a 3 = 0.35a 4 = 2.34(ft/i)- 1 /i- 1 Mpc (Bardeen et al 1986). On 
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scales greater than the horizon scale at the matter-radiation equality, ~ 13(£lh)~ 1 h~ 1 Mpc, the 
spectrum retains its original form and power slope index n. On smaller scales it is modified in a 
unique, for given Vth, way. The values of n that are usually considered are n ~ 1, required by the 
simplest CDM model based on inflation. However, a case has been made also for a tilted CDM 
model with n ~ 0.7 (Cen et al 1992). 

In models invoking topological defects, such as strings, the evolution of the initial density 
field can also be calculated although with significantly more uncertainty because of the possible 
evolutionary modes of strings when they enter horizon (Albrecht and Stebbins 1992). In 
models which involve only baryonic matter, the primeval density fluctuations have to be purely 
isocurvature, and the evolution of the density field is further complicated by the possible reheating 
and reionisation effects after recombination (e.g. Peebles 1987). 

In what follows we will distinguish between the "primordial spectrum" and what we term 
the "pregalactic spectrum" of the density field. By primordial we will mean the spectrum that 
was presumably produced during the very early stages of the Universe's evolution and which is 
assumed to be self-similar and characterised only by its power slope index n. "Pregalactic" in the 
terminology of the paper will refer to the spectrum of the density field after recombination but 
prior to galaxy formation, e.g. at redshifts > 100. In principle, the primordial spectrum can be 
recovered from the pregalactic one by assuming a transfer function. We will not do this, since it 
involves various uncertain assumptions, such as the nature of the dark matter, the type of density 
fluctuations, the values of cosmological parameters, assumptions about the degree and history of 
reionisation after recombination, etc. 

In what follows we will aim to reconstruct the pregalactic density field in the Universe from a 
variety of astronomical data, each dataset responsible for uncovering the spectrum of the density 
field over a certain range of scales. Because of the theoretical appeal of the CDM models we will 
compare the results in each range of scales to the anzatz eqs. (1),(2) predicted by inflationary 
cold-dark-matter models. 

In all models of the early Universe evolution, the density field on the largest scales which were 
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always outside the horizon during radiation dominated era, 3> a few hundred h~ Mpc, is expected 
to have preserved it initial, primordial, form. The spectrum of the density distribution on such 
scales is constrained by the COBE DMR data probing the microwave background anisotropies 
on angular scales > 7°. Assuming a scale- free primordial power spectrum with adiabatic initial 
conditions, COBE DMR data (Smoot et al 1992, Bennett et al 1996), after correcting for the 
Galactic cut, imply a roughly power-law spectrum with n ~ 1.1 ± 0.3 at the 68% confidence level 
(Gorski et al 1996). Normalization to the COBE data also fixes the value of the bias parameter 
b for a given power spectrum of the density field and f2 (Kashlinsky 1992; Efstathiou, Bond and 
White 1992). CDM models normalized to the COBE DMR data require b > 1 (e.g. Stompor et al 
1995). 

On smaller scales the pregalactic density field can be constrained and, as we show below, 
determined from the astronomical data out to z of a few. In the remainder of this article we 
recover the pregalactic density field from the various datasets on the present-day Universe. 

3. Two-point galaxy correlation function on linear scales: constraining the 
spectrum of light distribution on scales > 10/i _1 Mpc 

As was mentioned in the Introduction the two-point galaxy correlation function, £(r), and 

the power spectrum represent a pair of three-dimensional Fourier transforms and, for Gaussian 

density fluctuations, define all properties of the density field. They are related via: 

1 f'°° 

£(r) = ^2 I P(k)jo(kr)k 2 dk (3) 

Measurements of the two-point correlation function from galaxy catalogs show that on small scales 
£(r) is very close to a power law £(r) = (r/r*) -1-7 , where r* = 5.5/t _1 Mpc and 7 = 0.7 (Groth 
and Peebles 1977). This in turn means that the rms fluctuation in galaxy counts 

<(V^) 2 > = 3 / t(xr)x 2 dx = — 2 \ P(k)W TH (kr)k 2 dk (4) 
iVgai Jo 2ir z Jo 

is unity on scale r§ = (53^) ~ 8h Mpc (Davis and Peebles 1983). Here j n (x) is the 

n-th order spherical Bessel function and Wth = [3ji(x)/x] 2 . Eq. (4) also defines the connection 
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between the "counts-in-cells" analysis used in some galaxy surveys (e.g. Saunders et al 1991) and 
the underlying power spectrum or the 2-point correlation function, £(r). 

The data on the correlation function on larger scales are now available from the APM catalog 
measurements (Maddox et al 1990, hereafter MESL). The APM survey contains about 2.5 million 
galaxies in the blue magnitude range 17.5 < bj < 20.5. MESL measured the projected 2-point 
angular correlation, w(9), down to a systematic error of w pp ~ 1.5 x 10~ 3 , remaining due to 
plate-to-plate gradients. These data are consistent with a variety of measurements of £ from 
measurements at other bands (Picard 1991), catalogs (Collins, Nichol and Lumdsen 1992) and 
counts in cells analysis of IRAS galaxies (Saunders et al 1991). These measurements of w(0) probe 
£(r) to considerably larger scales than the non-linear scales of r < r§. 

If galaxy distribution along the line of sight is known one can relate w(9) to the 3-dimensional 
correlation function, £(r), via the Limber equation (Limber 1953; Peebles 1980). Its relativistic 
form is: 

w{6) = / dzftz) / dA£(r 12 ;z) (5) 

JO J-oc 

where r\ 2 = (c^|A) 2 + x ^+zf ^ s * ne P ro P er length, x(z) is the co-moving distance to z, 
cf)(z) = (dN I dz I N tot ) 2 accounts for the fraction of galaxies in the redshift interval [z;z + dz] and 
Ntot = Jo°° in~dz. The number of galaxies per dz with apparent magnitude mi < m < m u is given 
by: 

dN dV r L{ ~ m i) 

-^ = i~\ HL;z)dL (6) 
dz dz JUm,,) 



where dV/dz = ^+°^i+ nz f° r A = 0, &(L; z) is the galaxy luminosity function at redshift z and 
L(m) is the absolute luminosity of galaxy with apparent magnitude m at redshift z. 

Eq.(5) can be rewritten to relate w(9) directly to the power spectrum (Kashlinsky 1991a; 



Peacock 1991). Substituting (3) into (5) and using that Jf^ jo(\/x 2 + y 2 )dy = ttJq{x) leads to: 

/ °° dz(cdt/dz)-^(z) / °° dkPjk- z)Uq (iff) 
W[0) ~ 7T J °° P(k;0)W TH (kr 8 )k 2 dk (7) 

where Jo 0*0 is the zero-order cylindrical Bessel function. Note that deducing the 3-dimensional 

power spectrum from the projected galaxy angular correlation function is independent of the 
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distortions caused by peculiar velocities in the distribution of galaxies mapped in the redshift 
space (Kaiser 1987). Throughout the paper we will deal mostly with the shape of the power 
spectrum, rather than its amplitude. Hence we chose to normalize eq.(7) by introducing in the 
denominator the explicit expression, eq.(4), for unity galaxy counts over a sphere of radius r§. As 
eqs.(3),(5),(7) show a power law for w(9) oc #~ 7 implies a power law behaviour for £(r) oc r _1 ~ 7 
or power spectrum P{k) oc /c 7 ~ 2 . 

The APM survey represents the most accurate and biggest dataset for determining the 
correlation function galaxies. On small scales the projected angular correlation behaves like a 
power law w(9) oc 9~ 1 with slope 7 = 0.7. On larger scales w{9) falls off sharply and eventually 
becomes lost in the systematic errors dominated by the plate to plate gradients of the APM catalog, 
w pp ~ 1.5 x 10~ 3 . The falloff clearly implies the rollover in P{k) at large scales, or small k. The 
APM catalog contains galaxies down to the limiting blue magnitude of bj = 20.5 corresponding to 
z ~ 0.2. Thus in order to interpret the data on w{9) in terms of the power spectrum today (or at 
any other coeval epoch) one must eliminate or reduce the possible effects related to evolutionary 
uncertainties of the luminosity function, K-correction, and galaxy clustering. 

Various fits have been suggested for the power spectrum to describe the APM data. 
Kashlinsky (1992a; hereafter K92) suggested an empirical fit to the power spectrum to describe 
the data: on small scales P{k) oc k~<~ 2 to reproduce w{9) oc #~ 7 at small angles. On large scales, or 
k < ko, the power spectrum was assumed to go into the Harrison-Zeldovich regime consistent with 
the COBE DMR data. K92 used the APM data divided into six narrow magnitude of Am ~ 0.5. 
This way one can reduce effects of evolution since such narrow slices are more likely to contain 
galaxies at similar depths. Approximating the selection function like that of the Lick catalog 
(Groth and Peebles 1974; Peacock 1991) K92 found that satisfactory fits to the APM data can be 
obtained for k^ 1 ~ 40/i~ 1 Mpc. 

Baugh and Efstathiou (1993 - hereafter BE93) have developed an iterative deprojection 
technique using the Lucy method (Lucy 1974). They approximated the galaxy selection 
analytically and assumed that the time evolution of the power spectrum is scale independent. 
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BE93 then applied the method to invert w{0) from the entire APM dataset of 17.5 < bj < 20.5. 
The resultant three dimensional power spectrum with the error bars from the method and 
the APM data uncertainties is plotted in Fig. 7 of BE93. Application of their method to the 
two-dimensional power spectrum from the entire APM dataset led to similar numbers for P(k) 
(Baugh and Efstathiou 1994). 

In terms of CDM models, the APM data require significantly more large-scale power than 
predicted by the standard CDM model with Q = l,h = 0.5, n = 1 (MESL). Various suggestions 
have been made to account for the observed excess of large scale power on scales where the 
present-day density field is still linear, r > r 8 . Efstathiou et al (1990) have noted that the 
required large-scale power can be reproduced by low-0 CDM models in which case a non-vanishing 
cosmological constant would be required to keep the Universe flat in accordance with the standard 
inflationary scenario. Cen et al (1992) suggested that the excess power can be explained by 
the so-called tilted CDM model with n ~0. 7-0.8. Either of these modifications, decreasing VL or 
introducing n < 1, would boost the large-scale power of the CDM models on scales which are in 
the linear regime today. (At the same time this would, however, suppress the small scale power in 
such models; implications of this are discussed later in the paper). 

We now turn to quantifying how good the above fits and models for P(k) are viz-a-viz the 
various evolutionary corrections that can affect the accuracy of the various interpretations. The 
evolution of the luminosity function will affect the selection function of the Limber equation, 
while evolution of spectral energy distribution will determine which galaxies appear in the 
B-band at z = 0. Finally, the rate of the evolution of clustering pattern will introduce another 
redshift-dependent factor in the integrand of the numerator of (7). In addition, eq.(7) contains a 
dependence on O. 

MESL presented their data in two ways: 1) the entire sample of 2.5 million galaxies with bj 
magnitudes between 17.5 and 20.5, and 2) the data for sub-samples of galaxies binned in 6 narrow 
magnitude slices of Am ~ 0.5. In the first dataset, which was used in BE93, the galaxies are at 
very different depths and possibly very different evolutionary stages. The second APM dataset 
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used in K92 is more immune from effects of evolution since it is more likely to separate galaxies 
at different depths and any evolution can be more readily uncovered and corrected for. We thus 
consider the APM data on w(9) for the six slices each Am ~ 0.5 in width. 

Fig.l plots the selection function, zcj)(z), for the six slices with Am ~ 0.5 of the APM survey 
from bj = 17.5 to bj = 20.5. The luminosity function in eq.(6) was adopted from Loveday et al 
(1992): i.e. the Schechter (1976) luminosity function with M£ = —19.5 for h = 1 and a = —1. 
The relativistic K correction was modelled as 5m = Kz with K = 3 which adequately describes 
galactic spectra in the visible bands (cf. Yoshii and Takahara 1988). Solid lines correspond to 
SI = 0.1 and dotted to = 1 and no evolution was assumed. One can see that there is little 
dependence on Q. Similarly, there would be little change for other plausible values of < K < 4. 
The nearest slice of the APM galaxies contains galaxies typically lying at z ~ 0.07 and the farthest 
slice contains galaxies at z ~ 0.2, so certain amount of evolution could have occurred. There is 
no overlap at the FWHM level between galaxies in the nearest and most remote slices, although 
certain overlap exists between galaxies in the nearby slices. 

On co-moving scales less than r$ the present day density fluctuations are non-linear and the 
clustering pattern has been strongly distorted by gravitational evolution (e.g. Davis et al 1985). 
On larger scales the density field is still in the linear regime and the clustering pattern there is 
likely to have preserved the original power spectrum of the density field. We use the small scales, 
r < f8, where the power spectrum is well known and its evolution with z better understood, to 
constrain the possible evolutionary effects that can plague the APM data interpretation. 

As 9 — > the APM two-point correlation function is approximated well by a pure power 
law. At the same time, for small angular separations, the contribution to eq.(5) from £(r) at 
large separations becomes negligibly small. Thus in the limit of small angular separations the 
3-dimensional two-point correlation function at z = can be taken from the Lick catalog to 
be £(r) = (r/r*) -1-7 with r* = 5.5/i -1 Mpc. The time dependence of £(r) can be modelled in 
the very non-linear regime as £(r; z) = £(r; 0)^ 2 (z). Two extremes of clustering evolution are 
generally considered for very small scales (Peebles 1980): If clustering is stable in co-moving 
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coordinates ^> 2 (z) = (1 + z) 3 ; if it is stable in proper coordinates ^ 2 {z) = (1 + z) 7 . Thus for 
small scales the evolutionary effects must be constrained to match the APM data for all six slices: 
w{9) = A WtShcc 6-~t. Substituting £(r;z) = (r /r*)- 1 ^^^) into (5) leads to: 

7 

* 2 (z)(l + z) 2 y/l + Vtzdz (8) 

Comparison between (8) and the amplitude measured in the APM data provides an integral 
constraint on the amount of possible evolution. 

Fig.2 plots the data from the APM survey (MESL) divided into 6 narrow magnitude bins. 
The magnitude limits for each bin are shown on the top of each box. The data is plotted with 
open triangles. The power-law fits corresponding to AwO" 1 are plotted with two sets of straight 
lines. Dashed lines are for f2 = 1 and dotted are for O = 0.1 and A = 0. The two lines of each type 
correspond to clustering pattern stable in co- moving (lower lines) and proper coordinates. The 
amplitudes A w are shown for K = 3. One can see that no-evolution models describe the low-angle 
behaviour of w(9) extremely well. The dependence on can be neglected and the dependence on 
the precise rate of the evolution of clustering pattern is also small. For clarity of the figure we do 
not show the lines for other values of K, but note that the change is very small. Increasing the 
value of K would move the lines a little up, decreasing it to < 3 would make them lie exactly on 
top of the data at small angular scales. Similarly, making the clustering pattern evolve a little 
faster than ^ 2 {z) = (1 + z)~ 3 , as is suggested by simulations of non-linear gravitational clustering 
(Melott 1992) would shift the lines even closer to the data points. Thus we conclude that the 
APM data are consistent with no or little evolution out to the epochs accessible to the APM 
magnitude limit of bj = 20.5. The slight excess in the computed amplitude A w can be easily 
reduced to zero by adopting a slightly lower value of K or requiring the the clustering pattern to 
evolve faster than ty 2 (z) = (1 + z)~ 3 . The first of these is quite reasonable given the evolution 
of stellar populations in the relevant bands (e.g. Bruzual 1983). The second could be required if 
mergers play significant role in the evolution of clustering on non- linear scales (Melott 1992); this 
may even be suggested by the deep galaxy counts data (Broadhurst et al 1992). Thus we conclude 
that luminosity evolution plays a minor role out to z ~ 0.2 for the standard population of blue 
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galaxies with bj < 20.5. The no-evolution models with non-linear clustering stable in co-moving 
coordinates can be used sufficiently accurately in interpreting the APM data on w(9) in terms of 
the underlying power spectrum of galaxy distribution. 

Thin solid lines in Fig.2 show the BE93 fits to the six slices of the APM data assuming 
no luminosity evolution and with the luminosity function adopted from Loveday et al (1992). 
The three lines correspond to the best determined P (k) of BE93 and to one standard deviation 
uncertainty. The clustering pattern was assumed to be stable in co-moving coordinates on all 
scales. The BE93 power spectrum fits well the data at the various depths. It under-predicts the 
power on large scales by a small margin for the most distant of the APM slices. This, however, 
can be accounted for by assuming that the power spectrum on linear scales evolves less rapidly 
with redshift than the ty 2 (z) = (1 + z)~ 3 adopted from the non-linear scales. Indeed, linear scales 
should evolve differently in accordance with the growth predictions in linear regime. In principle, 
one can try to reproduce the consistent evolution of the power spectrum on all scales following 
the prescription of Hamilton et al (1991). However, it is not yet clear how well such fits work for 
arbitrary cosmologies and (non-power-law) power spectra (cf. Peacock and Dodds 1994). The fit 
of BE93 power spectrum to the APM data is good and we adopt it as a reasonable approximation 
to the power spectrum of galaxy distribution. The K92 fit is essentially the same as shown in 
Fig.l of the K92 paper and for brevity and clarity we do not show this set of lines in Fig.2, other 
than to point out that the BE93 power spectrum fits the APM data much better than the earlier 
K92 empirical fit. 

The K92 and BE93 power spectra fit the clustering hierarchy today over both linear and 
non-linear scales. On linear scales the power spectrum today should reflect initial conditions, but 
on scales less than rg the initial power spectrum was distorted by gravitational effects. CDM 
models predict power spectrum only in the linear regime so it makes little sense to plot them 
in the already condensed Fig.2. In order to constrain CDM models from the binned APM data 
we proceed as follows: the APM data is accurate for w(9) > w pp ~ 1.5 x 1CT 3 and the small 
angular scales analysis suggested that evolutionary effects are small for scales and depths probed 
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by the APM survey. The largest angular scales, where w{9) can still be probed, are dominated by 
contributions from linear scales, r > r 8 , and the data should be reproduced by eqs.(l),(2) within 
the framework of CDM models. 

Thus for a given CDM model specified by n and Qh we computed the angular scale, 9(w), on 
which the predicted angular correlation function drops to the value w. Fig.3 plots the resultant 
6(w) vs the mean magnitude, bj, of the slices in Fig. 2 for the APM data and for predictions of 
the standard (n = 1, upper boxes) and tilted (n = 0.7, lower boxes) CDM models. The left boxes 
are for w = 2 x 10~ 3 , which is just above the systematic error of w pp = 1.5 x 10~ 3 induced by 
the plate-to-plate gradients in the APM survey. The right boxes are for w = 5 x 10~ 3 which is 
significantly larger than w pp and should be quite accurate. The thick plus signs show the values of 
9(w) derived from the APM data. CDM models are shown with the other signs: x's correspond 
to Qh = 0.5, squares to Qh = 0.3, triangles to Qh = 0.2, rhombs by Qh = 0.1 and asterisks to 

Qh = 0.05. 

As was mentioned the values of 9(w = 5 x 10~ 3 ) are more reliable, but there is general 
agreement between the CDM fits to angular scales at both values of w. There is also consistency 
between the fits for various slices or depths. The fits in Fig.3 show that the linear part of the APM 
data can be described by CDM models which would then require Qh = 0.2 if n = 1 or Qh = 0.3 if 
n = 0.7. Increasing n above 0.7 would require even lower values of Qh. The limit on Qh required 
by CDM models for n = 1 is in good agreement with what was claimed before (Efstathiou et al 
1990). It is worth noting that the limits on the power parameter Qh for the CDM models come 
from the nearby slices. The differences in the values of 9(w = 5 x 10~ 3 ) between CDM models of 
different Qh are most profound for the nearest (and least affected by any evolution) slices. For the 
first three slices the APM data (pluses) would overlap with CDM predictions only for Qh = 0.2 if 
n = 1 and for Qh = 0.3 if n = 0.7; CDM models with both smaller and larger Qh would predict 
values of the angular scale 9(w) where w = 5x 10~ 3 ~ 3.3w pp different from what is observed. 
The limit on Qh (=0.3) for the tilted model is significantly lower than what was suggested before 
in Cen et al (1992) who used the entire APM galaxies lumped together and then scaled to the 
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depth of the Lick catalog. Thus even if the primordial n were as low as 0.7 one would still require 
low to fit the APM data with the CDM power spectrum. 

Strictly speaking, the above discussion applies only to the distribution of light (galaxies). 
In order to relate the power spectrum of galaxy distribution to that of mass one has to make 
assumptions about biasing or whether and how the light traces mass. The most economical and 
commonly made assumption to make is that of linear biasing (Kaiser 1984), i.e. that light traces 
mass at least to within a scale-independent constant. We followed this assumption in this section, 
but it must be remembered that the plots in Figs. 2, 3 refer to the distribution of luminous matter 
only. The next section discusses this assumption in more detail and there we propose an empirical 
justification for it based on comparing the density field of the APM survey to that probed by the 
velocity data. 

Assuming that light indeed traces mass, information on the initial (pregalactic) power 
spectrum in present-day galaxy catalogs is preserved only on scales where the density field is 
still in the linear regime, r > rs- Thus Fig. 4 plots the r.m.s. density fluctuation, A(M), over a 
sphere containing mass M vs the co-moving scale this mass subtends for the spectra required 
by the APM data. It is plotted in units on Ag, the r.m.s. fluctuation over a sphere of radius 
rs, and the numbers are shown for r > 10/i _1 Mpc where the galaxy density field is still in the 
linear regime. This ratio is used because in linear regime it is roughly epoch independent and also 
because for linear biasing the quantity A(M)/As recovered from the APM data is independent 
of the (constant) biasing factor b. Solid lines are the K92 power spectra with the transition 
scale k^ 1 = 35, 45, 55/i _1 Mpc from bottom up. Dashed lines correspond to the BE93 power 
spectrum with three lines showing the mean and ± one standard deviation uncertainty. The CDM 
models that fit the APM data in six narrow magnitude slices are shown with thick dashed lines: 
n = 1, Qh = 0.2 and n = 0.7, Clh = 0.3. (The two models would essentially overlap over the range 
of scales plotted in Fig.4). The various fits to the data practically overlap on scales < 50/i _1 Mpc. 
On larger scales the BE93 spectrum, which fits the APM data best, gives the least amount of 
large-scale power. Within the uncertainty of the methods the CDM models required by the APM 
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data reproduce roughly the same A(M)/A 8 ratio as the BE93 spectrum on scales < lOO/i x Mpc. 

The quantity A(M)/A§ plotted in Fig. 4 is independent of the bias factor provided the latter 
is scale-independent. It also measures the spectrum of the density field in linear regime and is 
approximately independent of the redshift at which it is evaluated. These are the main reasons why 
in this paper we choose to reconstruct this quantity from various datasets. 

To conclude this section: 1) By analysing APM catalog data in narrow magnitude slices it 
was demonstrated that evolutionary effects are unlikely to affect interpretation of the APM data. 
2) The power spectrum derived by BE93 provides a good fit to the APM data at all depths. 3) 
Assuming that light traces mass, the APM data can be fit at all depths by the standard CDM 
models with n = 1 if Qh ~ 0.2. 4) Tilted CDM models will require fi/i ~ 0.3 if n = 0.7 and lower 
values if n is increased. 5) In any case, all the above fits give practically identical numbers for 
A(M)/A 8 on linear scales < 50/i _1 Mpc. This quantity measures density field in the linear regime 
and is approximately redshift independent. If the possible bias factor is scale independent the 
quantity A(M)/As determined from the APM data is independent of the bias factor. 6) Fig. 4 
demonstrates consistency of this approach and we conclude that, provided that light-traces- mass, 
the pregalactic density field on these scales is close to that drawn in the figure. 

4. Peculiar velocity field vs galaxy correlation data: constraining VL and the bias 

factor 

The previous section established the power spectrum of galaxy distribution on scales 
10 < r <(60-100)/i _1 Mpc. In order to reliably identify it with the spectrum of the pregalactic 
density field one has to establish the relation between the light and mass distributions. 
Furthermore, as was discussed in the framework of the CDM models APM data would require 
a low O Universe. An independent test of these findings and assumptions can be provided 
by inter-comparison between the observed properties of peculiar velocity and density fields 
(Kashlinsky 1992b, 1994). We address these questions in the section. 
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Measurements of peculiar velocities provide one with direct probe of the mass distribution 
in the Universe and thus set strong limits on the large-scale structure models (e.g. Vittorio, 
Juszkiewicz and Davis 1987). Recent advances with the new distance indicators (Djorgovski and 
Davis 1987; Dressier et al 1987a) allowed one to measure the peculiar flows in the local part of 
the Universe up to ~ 50 — 100/i _1 Mpc (see Strauss and Willick 1995 for a review). The results 
from using the fundamental plane properties of elliptical galaxies (Dressier et al 1987b) suggest a 
large coherence length and amplitude of the local peculiar velocity field. The results from analysis 
of 1,355 spirals using the Tully-Fisher relation to determine distances (Mathewson et al 1992) 
would lead to an even larger coherence length and a similar amplitude. Other samples are also in 
agreement with the Great Attractor findings (e.g. Willick 1990). In a more controversial finding 
Postman and Lauer (1995) find from analysis of clusters of galaxies that the peculiar flows can 
be coherent over a much larger scale than the Great Attractor findings. All the results indicate 
significant deviations from the Hubble flow with a large coherence length. 

Since peculiar velocities probe the mass distribution and also depend on the density parameter 
£1 one can determine the latter by comparing the mass distribution implied by the velocity field 
with the observed distribution of galaxies. The density parameter is thus determined to within the 
uncertainty of the bias factor b by determining the factor Q,°- 6 /b. There are various ways to use 
this information to get at f2. Bertschinger and Dekel (1989) developed POTENT method whereby 
the mass distribution is reconstructed by using the analog of the Bernoulli equation for irrotational 
flows. They used the method to analyse in great detail and accuracy the peculiar velocity field out 
to about 60/i _1 Mpc (Bertschinger et al 1990). In Dekel et al (1993) they compared the previously 
determined velocity field with the observed distribution of galaxies and concluded that the best fit 
is achieved with Q os /b ~ 1. 

A different method was presented by Kashlinsky (1992b) and was based on comparing the 
velocity correlation function with the galaxy correlation function from the APM catalog; this led 
to low values of Q os /b ~0.2-0.3. Herebelow we present a more detailed analysis based on the latter 
scheme and show that the results are indeed suggestive of low values of the density parameter. We 
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also show that the results are in agreement with the scale independent bias factor. 

We define the dot velocity correlation function u{r) = (v(x + r) ■ v(x)) (Vittorio, Juszkiewicz 
and Davis 1987, Peebles 1987). 



Since gravity force is conservative, the flow caused by it must be irrotational. For irrotational flow 
the fc-th component of the Fourier transform of the peculiar velocity field is = —ikH f(Q)5k/k 2 
where f(£i) = dlnSk/dlnt. For the case of zero cosmological constant f(£i) = Q°' 6 /b. If the 
cosmological constant, A, is not zero f(Q) will be different, but for the most interesting case of 
U + A/3Hq = 1 the above approximation would work very well (Lahav et al 1991). Eqs (3), (9), (10) 
are also valid for the filtered fields in which case both v(r) and £(r) must be replaced with their 
filtered versions. 

Assuming that light traces mass to within a constant bias factor, b, substituting vu into (9), 
and then taking Laplacian operator of both sides would lead to: 



where it was used V 2 jo(A;r) = — /c 2 j (Ax) along with eq.(3). This equation was evaluated 
independently in Kashlinsky (1992b) and Juszkiewicz and Yahil (1989) and it can be used to relate 
the data on both peculiar velocity field and galaxy correlation function to get f2°' 6 /& (Kashlinsky 
1992b, 1994). The advantage of using (10) rather than the power spectrum data in conjunction 
with (9) (e.g. Kaiser 1983) is that eq.(10) can be integrated to contain only the scales over which 
the data on both £(r) and v(r) are known, while in the former case assumptions have to be made 
about the power spectrum over the wavelengths inaccessible to observations (Vittorio et al 1987). 
(Of course, the data £(r) are in effect on integral constraint on P(k) over the infinite range of 
wavelengths). Equation (10) can be solved to give the velocity correlation function in terms of 
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v{r) = i/(0) - 
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where it is defined J n (r) = /J* ^(x)x n ~ 1 dx. Only one of the constants of integration remains in 
(11), the other constant of integration vanishes because the velocity field must be finite at r = 0. 
The constant of integration remaining in (11) is just the "central" velocity dispersion, u(0) = (v 2 ). 

Strictly speaking in order to use eq.(10) the data on both u(r) and £(r) must come from the 
same region of the Universe. Since the current velocity data sample galaxies within a radius of 
~ 100/i~ 1 Mpc, the numbers for £(r) must come from the "local" surveys such as the CfA survey 
(Geller and Huchra 1989). The CfA redshift survey is complete out to tub < 15.5 and for the 
relevant scales, r ~ 40ft,~ 1 Mpc, the correlation function in the CfA survey essentially coincides 
with the one determined from the APM catalog (Da Costa et al 1994; Vogeley et al 1992). In fact, 
the numbers in this section are in good agreement with each other when evaluated for £(r) from 
the CfA survey, the power-law model, or the BE93 or K92 fits to the APM data on w{9). The 
galaxy correlation function on the relevant scales from the APM data, as discussed in the previous 
section, is plotted in Fig.ll; for the CfA survey £(r) is plotted in Fig.3 of Vogeley et al (1992). 

The observed correlation function of galaxies would predict large velocities if J7 6 /6 = 1. For 
power-law f(r) = (r/r*)- 1 "^ eq.(ll) would lead to 0/(0) - v{r) = H r» n°- 6 (^§(1-7) = 

V( 1 -7)(2-7) * 

881^(^) ai5 km/sec for 7 = 0.7. The value of y/v(0) - v(r) for n os /b = 1 would be significantly 
larger than the data on peculiar velocities at 40 or 60/i _1 Mpc of Bertschinger et al (1990) (see 
Table 1). Fig.5 plots the predicted value of y/ v(0) — v{r) computed according to eq.(ll) for 
Q°- 6 /b = 1; the numbers in Fig.5 scale oc Cl°' e /b. The data on £(r) from which the lines in Fig.5 
were computed are for the APM correlation function using the BE93 deprojection (shown with 
solid lines) and the correlation function of the CfA survey from Vogeley et al (1992) (dashed line). 
The three solid lines correspond to one standard deviation uncertainty in BE93. CfA survey 
measures £(r) out to only ~ 40ft,~ 1 Mpc; hence the line ends there. The lines show that there is 
good agreement with the dot velocity correlation function values computed "locally" (CfA) and 
"globally" (APM). There is also good agreement between the shape of the bulk velocity flows 
in the Great Attractor region (Bertschinger et al 1990) and Fig.5. This is consistent with the 
assumption of scale independent bias parameter b. 
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We ask now the following question: given the measurement of v(r) on some scale what should 
one expect to find for i/(0) given the data on £(r)? This can be evaluated from 

1/(0) = i/(r) + ^flg[J 2 (r) - ^] (12) 

In order to evaluate z/(0) from i/(r) we use the numbers from the POTENT analysis (Bertschinger 
et al 1990) who compute the bulk velocity V(r) field after filtering the data with a Gaussian filter 
of filtering length rj = 12/i _1 Mpc. The numbers from their analysis of the velocity field are given 
in Table 1 which shows the values of V(r), the amplitude and direction, for three different r. 

The data with which i/(0) ought to be compared are admittedly incomplete, but it is generally 
thought that for non-filtered velocity field y/vQS) ~ (500-600)km/sec (cf. Peebles 1987). This 
also follows from the numbers for velocity dispersions in typical collapsed systems, such as groups 
of galaxies. Local measure of y/v{fy is probably the dipole velocity which is observed to be 630 
km/sec with a very small error. (The direction of the microwave background dipole roughly 
coincides with the velocity vector from Bertschinger et al (1990) in Table 1.). Since v(r) — > as 
r — > oo, eq.(12) shows that the linear perturbative expression for u(0) is oc J2{r — > oo) = J2,oo- 
In this case, because £(r) oc r -1 ' 7 at small scales, a non-negligible contribution to z/(0) comes 
from non-linear scales. This suggests the importance of smoothing when evaluating v(0) directly 
from the data. Peebles (1988) utilized the exact (non-perturbative) Layzer-Irvine cosmic energy 
equation to compute z/(0) using the value of J2 )0 o = 164e ±ai5 /i~ 2 Mpc 2 evaluated from the Lick 
data (Clutton-Brock and Peebles 1981). He obtained y/u(0) ~(500-600) km/sec for in the 
range 0.2-0.3 in reasonable agreement with this discussion. For comparison, integrating over £(r) 
from the BE93 deprojection of the APM data gives </2,oc ^ 170 ± 5h~ 2 Mpc 2 , while for the CfA 
data from Vogeley et al (1992) one obtains J2,oo — 186/i~ 2 Mpc 2 . Both are in agreement with the 
numbers used in Peebles' (1988) analysis. 

Table 2 shows the values of \J u(0) evaluated from eq.(12) given the velocity data input on 
v(r) ~ V 2 (r) at r =40 and 60/i~ 1 Mpc from Bertschinger et al (1990). The columns show numbers 
for £(r) on scales < 60/i _1 Mpc computed according to pure power- law, £(r) = (r/r*) -1-7 , the K92 
and BE93 fits to the APM data; and the CfA data on £(r) from Vogeley (1992) at r < 40/i~ 1 Mpc. 
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The first column gives the filtering radius in h~ 1 Mpc, the second the value of Q os /b used in (3), 
the next seven columns give the values of \J u(0) in km/sec computed according to eq.(12) for the 
various The first four rows show the numbers with no filtering, the last four rows show the 

numbers for rj = 12/t _1 Mpc used in Bertschinger et al (1990) analysis of the velocity field. The 
data for the CfA measured £(r) are shown only for r = 40/i~ 1 Mpc, the scale up to which the CfA 
data can probe £(r). We have not shown the filtered values for the CfA £(r) because the latter 
was not measured out to large enough scales to make filtering with r/ = 12/i _1 Mpc meaningful. 

The various fits to the galaxy correlation give similar numbers which shows the robustness of 
the various methods. The data on £(r) determined from the CfA survey come from the region 
that includes galaxies used in the peculiar velocity analyses (Faber et al 1989) and agrees well 
with z/(0) predicted by the APM data. This is turn means that the Great Attractor is a typical, 
rather than rare, mass concentration in the Universe. There is very good agreement between 
the numbers in Table 2 evaluated at both on 40 and on 60/i _1 Mpc. This is consistent with the 
constant bias factor over at least this range of scales. 

Low values of Q are implied by this analysis. E.g. for the non-filtered field the predicted 
value of y/v(0) should be compared with the data on the dipole velocity of only 600 km/sec. Also 
if S7°- 6 /6 were as high as unity the typical velocity dispersion in the collapsed systems would be 
around 1300km/sec, corresponding to X-ray temperature of >10 Kev. Observations, however, 
suggest that typical collapsed structures, such as groups and poor clusters of galaxies, have 
velocity dispersion ~ 500 km/sec and rich X-ray emitting clusters contain only a small fraction 
all galaxies. For the filtered field with rj = 12/i _1 Mpc the predictions should be compared with 
the value from the Bertschinger et al (1990) analysis that gives 457 ± 61 km/sec. Again the best 
agreement between the data and u(0) would be achieved for O 6 /6 ~ 0.2 — 0.3; the fact that we 
reach essentially the same conclusions when comparing the filtered v{r) is reassuring. 

Thus the above analysis suggests that low-fi Universe may not be in conflict with the data 
on peculiar velocities and may in fact be even suggested by the latter. Similar conclusions are 
reached applying the least action method (Peebles 1989,1990) to the dynamics of the local Group 
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(Shaya et al 1995). We have not quantified the conclusions from Table 2 in statistical terms and, 
hence, take them only as suggestive. Below we propose a further extension of this discussion that 
may determine the values of 17 6 /6 in the new datasets on peculiar flows. 

A similar analysis can be applied to the components of the velocity tensor: 
Uij = (vi(x + r) ■ Vj(x)) = T>(r)5ij + [n(r) — S(r)] '.' 3 . For irrotational fluid the transverse 

% 3 

velocity correlation function, £(r), and the parallel, Il(r) are interrelated via: 

3rX(r) 



n(r) 



dr 



and the dot (total) velocity correlation is given by (e.g. Gorski 1988) 



(13) 



n(r) +2S(r) = v{r) 



(14) 



Eqs(10),(13) and(14) allow one to evaluate both components: 
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Now note that the quantity £(r) — II(r) is independent of the (a priori unknown) integration 
constant v(0) and hence, once measured, it can provide a tool for measuring O 6 /6. It is given by: 

,^ L2 flJ 3 (r) U 5 (r)- 



E(r) - n(r) = Hl-p 



2 r 3 



(17) 



Since the logarithmic slope of £(r) is > —2, the expression on the right-hand-side of (17) is 
dominated by the contribution from the integrands in J3, J5 near r. Therefore for this quantity it 
is sufficient to choose a large enough separation to ensure the validity of the linear approximation. 

For a power-law £(r) = (r/r*)- 1 ^ one gets ^S(r) - U{r) = ^ ^ (2 f° } r * 4 _ 7) (^)^ (1 ~ 7) ^(400- 
500) km/sec for 7 = 0.7 if S7 6 /6 ~ 1. Fig. 6 shows the predicted \/S(r) — II(r) in km/sec for 
0°' 6 /6 = 1. The lines were evaluated using the CfA and APM data on £(r); the notation is the 
same as in Fig. 5. Again there is good agreement between the values computed from the "locally" 
determined correlation function (CfA - dashed line) and the "global" £(r) from the APM data. 
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The amplitude in Fig.6 scales oc S7°- 6 /6 and is quite significant if Q os /b ~ 1. The data from 
previous analyses trying to reconstruct the velocity correlation tensor (Gorski et al 1989; Groth 
et al 1989) from, by now old, catalogs are too uncertain to use in comparison with Fig.6. But 
the magnitude of ^S(r) - n(r) of ~400 km/sec over the range of (30-60) h _1 Mpc if tt 0M /b ~ 1 
suggests that this could be a measurable task with the new datasets and that in conjunction with 
the galaxy correlation data this may further constrain Q 0M /b. 

Thus the analysis and results of this section suggest the following: 1) the spectrum of mass 
distribution from velocity data is consistent with that determined from the galaxy correlation 
data (cf. Kashlinsky 1992b). 2) This in turn suggests that the bias factor is consistent with being 
scale-independent. 3) Low values of O a6 /6 may be consistent with the data on peculiar velocities. 
4) Thus we interpret the results in Fig. 4 as the rms density fluctuation of the pregalactic mass 
density field on scales rs < r < 100/i _1 Mpc. 

5. Cluster correlation function and its richness dependence: reconstructing 
pregalactic spectrum on scales 5-25/i x Mpc 

The cluster-cluster correlation function is known to have a larger amplitude than the one 
measured by galaxies and its amplitude also increases with the cluster richness (Bahcall and 
Soneira 1983). Kaiser (1984) suggested that the increase in the amplitude can be explained 
if clusters formed at the rare peaks of the initial density field. His analysis was later applied 
by Bardeen et al (1986) to study in great detail the properties of peaks in the density field. 
Kashlinsky (1987,1991b) proposed that one can combine Kaiser's original suggestion with 
gravitational clustering theory (Press and Schechter 1974) to explain both the amplification and its 
correlation with the cluster richness. He supposed that structure formation in the Universe forms 
via a "natural bias" (cf. Davis et al 1985), i.e. systems of galaxies (clusters and groups) should 
be identified with regions that had turn-around time time less than the age of the Universe. The 
model successfully explained the correlation between the cluster correlation function and richness 
(mass) with objects that turned-around on a larger scale having a greater correlation amplitude. 
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Within the gravitational clustering model the properties of the hierarchy would then depend 
uniquely upon the initial power spectrum and the mass of the objects that formed/turned-around. 
In this section we first evaluate in greater detail the predicted properties of the cluster correlation 
function according to the gravitational clustering picture (Kashlinsky 1987, 1991b) and then apply 
the results to invert the data on cluster correlation function to obtain the pregalactic spectrum of 
the density field on scales 5/i _1 Mpc < r < 20/i~ 1 Mpc. 

We start with the density field at some initial epoch, S> 1, when density fluctuations 
were linear on all scales of interest. The final results will be independent of Z{. We define with 
8m = J 8(x)dV/ J dV the initial mass over-density over the comoving volume V that contains 
mass M. On large scales, r 3> R(M), where R{M) is the comoving scale containing mass M, the 
correlation function of the 5m field coincides with £(r). At zerodag it equals the mean square 
fluctuation over mass M: Af(M) = (5 M ); the subscript i refers to the values evaluated at Zj. Thus 
the correlation matrix for the <5M-field is: 



{Clm} 



Af(M) &(r) 
v fc(r) Af(M) 



(18) 



We assume that the initial density field was gaussian. Then the probability density to find 
two regions containing masses M\,Mi that had density fluctuations <$i,<$2 at Zi is: 

p{Sl ' S2) = 2^C\\ eM ~\ 6 ' C ~ K6) (19) 
As the clustering evolves deviations from gaussianity will develop due to gravitational effects. 
However, since we trace the distribution of the prospective clusters at z\ these will not be 
important for the computation of the cluster correlation function below. Using properties of the 
Fourier transform of a multidimensional Gaussian, eq.(19) can be rewritten in terms of the direct 
matrix C as: 

1 rco roc 1 

p(S 1 ;S 2 ) = ^ J J exp(-iq ■ 5)exp(--g • C ■ q)d 2 q (20) 

We denote with 5t a the amplitude at z% required for the fluctuation to turn around at z = 
and use the Press-Schechter (1974) prescription for gravitational clustering which assumes that 
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any region with 5m > $ta would have turned around by now. We further assume that clusters 
and groups of galaxies are identified with such regions. Then the probability for two clusters of 
masses M±, M 2 to have formed at any time between now and Zi is given by integrating (2) over all 
fluctuations > St a - We denote with Ag,i the amplitude the fluctuation had to have on scale r§ at 
Zi in order to grow to the observed value of 1/6 at z = 0. A convenient way to characterize 5t a 
and to normalize the density field is by introducing a quantity Qt a = ^f" (Kashlinsky 1991b). 
For A 8 = 1 today, Q ta — 1.65 almost independently of Q or Zj. We adopt this value of Q ta in the 
calculations in this section. We will show that this set of "minimal assumptions", including 6=1, 
is justified since it gives a fit to the derived spectrum consistent with the APM data. 



Expanding in eq.(20) exp( — ^q • C • q) = J2r 
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we can write the probability for two clusters of mass Mi , M 2 to form at any time between Zi and 
z = as: 



Pm 1 m 2 
where: 
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The quantity C is the number of standard deviations the mass M is with respect to the "typical" 
member of the hierarchy and is given by: 



C = Q 
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(23) 



/ °° P(k)W TH (kR(M))k 2 dk 

where R(M) is the scale containing mass M. The power spectrum in (23) is the original pregalactic 
power spectrum of the density field. 

In order to directly evaluate Pm\M 2 an d the cluster correlation function we write 
I-oc exp(-ixy - \y 2 )y m dy = i m -§^r JZo exp(-ixy - \y 2 )dy = i m V2^-§^ exp(-±x 2 ) (Jensen and 
Szalay 1986). Then eqs.(21),(22) become: 
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The probability of cluster of mass M to form by now is for Gaussian ensemble given by: 



(25) 



Eqs.(24),(25) give the cumulative probability of objects to turn around at any time in the 
past. In the gravitational clustering picture -Pm's include objects that formed at earlier epochs and 
by now are incorporated into larger systems. Because clusters and groups form dissipationlessly, 
the objects that formed earlier would merge into larger systems and lose their identity as clustering 
progresses to larger masses. Since we observe at a fixed epoch (z = 0), one must translate (24), (25) 
into probabilities that object formed (turned-around) today. Thus the fraction of pairs of clusters 
that formed today on mass scales M\ , M2 separated by distance r out of ensemble of density 
fluctuations given by eq.(19) is Jm 1 m 2 = qmom 2 • Similarly the fraction of clusters in the mass 



range [M,M + dM] is f M dM = j^fdM (Press and Schechter 1974). From (24), (25) one can show 
that 



where B m (x) = ^x m_1 H m+ i(x) (Kashlinsky 1991b) and H n (x) = (— ) n exp(x 2 )J^r exp(— x 2 ) are 
the Hermite polynomials. We will be using the measurements of the cluster correlation function 
on linear scales, r > r$. On these scales the ratio on the right-hand-side of (26), £j/Ag, where 
both the numerator and denominator are evaluated at Zi, is roughly independent of redshift and 
is approximately equal to the galaxy correlation function today, £(r). 

Having fixed the fractions of objects that turned-around today on a given scale we can 
evaluate the correlation function between clusters of different masses £miM 2 - By definition, 
the 2-point correlation function of an ensemble of objects with number density n is given 
by the probability to find two objects in small volumes dVi,dV2 as dV\2 = n 2 (l + £)dV\dV2- 
Since the clusters of mass M\,M.2 will make the fraction /miM 2 of such pairs, the probability 
to find them is dVM±M 2 = fM\M 2 dV\2- On the other hand, the mean number density of 
clusters of mass M would be fu x n and by definition the probability to find two clusters is 
dT > M 1 M 2 = /Mi/M2 ri2 (l + £,M 1 M 2 )dV\dV2- Hence the correlation function of clusters of mass M is 




(26) 



-29- 



given by: 



H IVJ 2 



1 + imm = iSt 2 {1 + ^ = vMk {1 + 6) (27) 



Eq.(27) fixes the factor by which the cluster correlation function is amplified over the 
underlying correlation function of the hierarchy, £. The amplification is purely statistical and the 
discussion does not involve any dynamical effects. Hence, on linear scales today the amplification 
factor £mm/£« is redshift independent and for r > r§ the present-day correlation function for 
clusters of mass M should be amplified over £(r) according to: 

= A M {^=)i(r) (28) 



where: 



A «W = E tfS^FCMix) (29) 

m=0 ^ta 
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(30) 



and 

x 2m "j_r2 

Cm ^ = — [- X 2 ' ( TO + 1)<& 
In evaluating (28)- (30) we used that the ratio ^/Ag^ evaluated at Zi for scales r > r§ is ~ 
the galaxy 2-point correlation function that is measured today. As x — > eq.(30) gives Co — > 1 and 
one gets from the first term in (28) that at small x: ^mmH — £(?")■ F° r large masses or s > 1 the 
first term in (29) is Co = ( 4 and the amplification reduces to ^mmH — Qta[^-s/^-(M)] 4 ^(r) (cf. 
Kashlinsky 1987). (The quantity Q ta in this discussion is equivalent to the threshold amplitude b 
of Kaiser 1984). 

Eqs.(28)-(30) show that the spectrum of the pregalactic density field can be constrained, and 
as we show later determined, by the data on both the slope of the observed cluster correlation 
function £cc(r) via eq.(28), and the dependence of the amplitude on the cluster mass via eq.(29). 
Bahcall and Soneira (1983) present the data on both of these for scales > r§ where the present 
analysis applies. The slope of the cluster correlation function on very large scales, r > 50/i _1 Mpc, 
is poorly determined from the data (cf. Fig.9 of Bahcall and Soneira 1983). For richness class 
1Z > 1 Abell clusters, such as Coma, they approximate £cc( r ) — 300(r/l/i _1 Mpc) -1 - 8 . This 
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approximation has large uncertainty at large scales, but for r < (30-40) /i _1 Mpc it can be used as 
a reasonable approximation to the data. The Zwicky clusters, which are poorer, also exhibit a 
stronger correlation amplitude than galaxies (Postman et al 1986), although it is weaker than that 
of Abell clusters and is consistent with the amplitude-richness relation proposed in Bahcall and 
Soneira (1983). Fig.2 in Bahcall and West (1992) shows the largest compilation of the data on 
cluster correlation amplitude vs richness. Thus we use below the data on the cluster correlation 
function from Bahcall and Soneira (1983) and on the amplitude-correlation amplitude dependence 
for clusters of galaxies from Bahcall and West (1992). 

In order to compare eqs.(28)-(30) to the data discussed above, we must fix the mass of the 
clusters quantitatively. We use Coma, the best studied cluster, as the mass normalization point. 
From Kent and Gunn (1981) we adopt its mass to be -/^fcoma — 

1.45 x 1O 15 ^, _1 M (see also White 
et al 1993 and references cited therein). The richness of the Coma cluster is adopted from Bahcall 
(1981) and Abell et al (1989) to be A/boma = 106. Assuming that both the luminosity function of 
galaxies and the ratio of dark-to-luminous matter are universal for all clusters would imply that 
the cluster richness is proportional to mass. We thus adopt the following relation between cluster 
mass and richness: 



while the comoving scale containing mass M that enters in the integrand of denominator of (23) 
is given by: 



We can now move to the implications of the cluster correlation data for the spectrum of the 
pregalactic density field. For a given P(k) eqs.(3) and (23) specify £(r) and A(M)/A§ which in 
turn uniquely determine both the cluster correlation function vs r and its amplitude vs the mass 
computed from the cluster richness according to eq.(31). For CDM models the value of Qh is not 
the only parameter that fixes the cluster correlation function; the extra dependence on h comes 
from translating the mass-scales to linear scales via eq.(32). Fig.8 plots the cluster correlation 
in CDM models with the Harrison-Zeldovich power spectrum, n = 1 for various Q and h. Thin 




(31) 




(32) 



-31 - 



solid lines correspond to the underlying correlation function, £(r), for given Q,h. Thick solid lines 
correspond to the £cc = 300(r/lft,^ 1 Mpc)~ 1 - 8 that Bahcall and Soneira find for TZ > 1 clusters. 
The numbers are plotted between 20/i _1 and 100/t _1 Mpc where the galaxy correlation function is 
in the linear regime and thus is given by the pregalactic power spectrum. Dotted lines correspond 
to Af = 5, dashes to M = 10, dashed-dotted lines to Af = 50 and dashed-dotted-dotted lines to 
clusters such as Coma, Af = 100. Successful models should fit both the amplitude and the slope 
of the measured £,cCi while at the same time the underlying £(r) should be in agreement with 
the APM data constraints in Fig. 3. As one can see from the figure, the required increase in the 
amplitude for clusters of richness class 1Z > 1 (Af ~ 100) can be achieved for ~ (0.3-0.4), but 
such models would not reproduce the observed slope of £cc(r) even on scales < (30-40) /i _1 Mpc 
where the data is more reliable and is definitely not expected to be subject to possible projection 
effects (Sutherland 1988). The required slope of icc( r ) can be reproduced for f2 = 0.1 and h = 0.5, 
but then the observed amplitude would be reached for Af = 10, corresponding to poor groups, 
while clusters such as Coma should have correlation amplitude significantly above that observed. 
Furthermore, such models with Qh ~ 0.05 would have values for 6(w = 5 x 10~ 3 ) - plotted with 
asterisks in Fig.3 - that are much bigger than the values deduced from the APM data. For brevity 
we do not present the same graph for h = 1, but the agreement with the measurements would 
become even worse for h > 0.5. 

Fig.8 plots the predicted ^mmW for tilted CDM models with n = 0.7; the line notation is the 
same as in Fig. 7. The tilted CDM model also does not fit well the data on the cluster correlation 
function or the APM data. E.g. the model can reasonably reproduce the slope and the amplitude 
of Abell clusters if Q = 0.3, h = 0.5, but then with Qh ~ 0.15 it would predict 9(w = 5 x 10 -3 ) in 
Fig.3 far above the APM data. 

In order to illustrate the dependence of the amplitude on mass or richness for CDM models 
we computed the amplification factor, eq.(29), at r = 25ftT 1 Mpc. This scale is » r§ so it is 
reasonable to suppose that the density field there reflects the pregalactic density field. At the same 
time, on this scale the underlying £(r) is still positive in all relevant CDM models. Fig.9 plots the 
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values of A(25h~ 1 Mpc) vs the cluster richness for CDM models. Right box shows the numbers 
for n = 1 and the left box shows them for tilted CDM models. Triangles correspond to the data 
from Fig. 3 of Bahcall and West (1992). Solid lines correspond to £lh = 0.1, dotted to flh = 0.2, 
dashes to flh = 0.3 and the dashed-dotted lines correspond to flh = 0.5. Thin lines of each type 
are for h = 1 and thick lines are for h = 0.5. One can see that it is difficult to describe the data 
with one particular CDM model. In other words, the slope of the pregalactic power spectrum over 
different range of scales as probed by Figs. 7,8, 9 cannot be fitted with one formula given by eq.(2) 
for any value of f2, h or n. Further difficulty would come from constraining the CDM models to 
fit the APM data along with the cluster correlation data. As Fig. 3 shows, CDM models would be 
consistent with APM data only if VLh ~ 0.2 for n = 1 and Qh ~ 0.3 if n = 0.7. This would further 
restrict the models to dotted lines in the left box of the figure and to dashed lines in the right box. 

One can now reverse the problem and instead of fitting various theoretical models to the 
cluster correlation amplitude - richness data, one can invert the data to obtain the implied 
pregalactic density spectrum from eq.(29) with x given by eq.(23). This can be done if one uses 
independent measurement of the underlying correlation function, £(r), in linear regime. Then the 
latter can be substituted into (29) with x given by (23) to give an equation solving which one 
can determine, given the measurements of A(M) vs richness for clusters, the values of A(M)/Ag 
of the pregalactic density field on mass scale M (in turn given by eq.31). First to the choice of 
scale r on which the data on £(r) is to be used. Such scale must be > r§ since £(r) in eq.(29) 
is given by the pregalactic power spectrum at z,. Fig.10 plots £(r) from the various fits to the 
APM data. Solid lines are the K92 fits with k 1 = 35, 40, 50/t _1 Mpc from bottom to top. Dotted 
lines are for the BE93 fit with top and bottom lines corresponding to the one standard deviation 
uncertainty. For comparison we also plot with dashed lines the correlation function according to 
CDM models (eq.2) with n = 1 for Qh = 0.2. (n = 0.7 and flh = 0.3 line would essentially coincide 
with the thick line over the linear scales; hence it is not plotted). All the fits to the APM data 
give essentially the same £(r) between 10 and 30ft, -1 Mpc. Thus we choose scale r = 25/i _1 Mpc to 
substitute the measured £(r) into (29). The latter is sufficiently larger than r§ to be certain that 
£(r) there reflects the pregalactic density field and, at the same time, sufficiently small so that 
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the cluster correlation amplitude measurements are free from the possible uncertainties at larger 
scales. Using the data on r anywhere between 15 and 30/i -1 Mpc would introduce little difference 
in the results to come. From Fig. 10 we adopt scale r = 25/i _1 Mpc with the measured correlation 
function there £(25/i _1 Mpc) = 0.07. The same number would be given also by a simple power 
extrapolation of £(r) = (r/r*) -1 - 7 to r = 25/i~ 1 Mpc. 

Solid line in Fig. 11 plots the amplification factor A(x) at 25/i _1 Mpc vs x according to 
eq.(29) after adopting £(25/t _1 Mpc) = 0.07. One can see that at small x (small mass scales) 
A(x) has a very weak dependence on x. Hence for small mass scales the determination of x, or 
A(M)/A 8 = is less accurate since the errors in the data will amplify. On the other hand, 
at large x the dependence is quite steep and the determined spectrum will be less sensitive to 
observational errors in cluster correlation measurements for massive clusters. Note that the first 
(and leading at large x and small £) term in (29) is independent of £(r). The two dotted lines 
show the uncertainty in A{x) introduced by assuming e.g. a 25 % uncertainty in the value of £, 
i.e. the lines plotted cover the range of 0.05 < £(25/i -1 Mpc) < 0.1. One can see that the values of 
A(M)/As determined in this way depend very weakly on the possible uncertainties in £(r). 

We now use the data from Fig.3 of Bahcall and West (1992), plotted with triangles in Fig.9, 
in order to invert eq.(29), with £(25/t _1 Mpc)=0.07, to give A(M)/As of the pregalactic density 
field on scale M, or richness M. Fig. 12 plots the results of this inversion. The top horizonal 
axis plots the values of jV at which A(M)/As has been evaluated. The bottom horizontal axis 
shows the mass computed according to the normalization given by (31). We emphasize again 
that this method gives directly the pregalactic spectrum as it was at Zi independently of the later 
gravitational or other effects. The plot in Fig.12 shows a clearly defined slope of A(M) oc M~ - 275 
corresponding to the spectral index of n ~ —1.3. This slope is consistent with the APM implied 
power spectrum index of w{6) oc 9~ - 7 . 

Finally, we note that even if the underlying density field is purely gaussian, the distribution of 
clusters of galaxies would be non-gaussian (Politzer and Wise 1984, Kashlinsky 1991b). Kashlinsky 
(1991b) analysed the properties of the three-point correlation function of clusters predicted by 



-34- 



this model and found that they compare favorably with the measurements of the cluster 3-point 
correlation function by Toth et al (1989) assuming that the pregalactic power spectrum is close 
to n ~ — 1 on all scales < 100/i _1 Mpc. In principle, one could apply similar methods to derive 
the properties of the pregalactic density field from the 3-point correlation function of clusters of 
various richness/mass. The data available at present do not justify a lengthy addition on this to 
the paper and we postpone this part of discussion to a forthcoming paper. 

6. Combining the results from w(9) and £ cc : Q, and the pregalactic density field 

The agreement of the slope of the spectrum in Fig. 12 determined from the cluster correlation 
amplitude vs richness data with that deduced from the APM catalog is encouraging and shows 
consistency of this approach. Furthermore, the APM dataset fixes the spectrum of the pregalactic 
density on a given linear scale, r, while the method outlined in the previous section determines it 
on a given mass scale, M. The resultant spectra are consistent in slope and requiring them to give 
identical amplitude at a given linear (or mass) scale could then constrain fi according to eq.(32). 

Fig. 13 shows A(M)/A$ from Fig. 12 plotted as function of the linear scale containing the mass 
on the lower horizontal axis of Fig.12. The linear scale r in h~ 1 Mpc was computed according to 
eq.(32). Left box of Fig. 13 shows the numbers for Q = 1 (thick plus signs) and Q = 0.1 (thin plus 
signs). The open square in the figure shows the value of 1 for A(M)/A§ at r§; for the pregalactic 
density field this number must be unity by definition. Clearly = 1 would be difficult to reconcile 
with A(M)/Ag = 1 at r$; extrapolating from all the data points misses the square by a significant 
factor. A value of Q < 1 would be required to reproduce the unity value for A(M)/A§ at r§ as 
obtained from the cluster correlation amplitude vs richness dependence. 

We now use (32) to determine O by requiring the pregalactic density field derived from the 
cluster correlation function to match that derived from the APM data on the same scales. The 
right box of Fig. 13 shows the intercomparison. The lines are in the same notation as in Fig. 10: 
solid lines correspond to the K92 fit and dotted lines correspond to the BE93 fit. The plus and 
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square symbols correspond to the same notation as in the left box of the figure. The rhombs 
correspond to 17 = 0.25 with which the best agreement between the A(M)/A 8 values derived from 
the cluster correlation amplitude and APM dataset is achieved. Note that any deviations from 
spherical approximation will speed up the collapse (Peebles 1980) leading to smaller value of Qt a - 
This will shift the values of A(M)/Ag from the cluster correlation analysis that fit APM data best 
towards even lower values of fi. 

It is interesting to note that the value of Q = 0.25 deduced in this analysis is in good 
agreement with that required by the dynamics of the Coma cluster. This cluster is by far the best 
studied and was also used in the normalization of the mass-richness relation used in Sec. 5. The 
blue mass-to-light ratio of the Coma cluster is very accurately determined to be M/Lb — 362/j in 
agreement with the value of Mcoma used in (31). If one adopts the blue luminosity function of the 
Schechter form $(L)dL = §*(L/L*)~ a exp(—L/L*)d(L/L*) and assumes that the entire mass in 
the Universe is associated with galaxies and their systems, the total mass-to-light ratio of galactic 
systems would be uniquely related to Q via: 

— = — ^ -n = 1330/i n solar units (33) 

L B 87rG$*L*r(2 - a) v ; 

where T(x) is the gamma function. In eq.(33) the numbers were evaluated using the data on the 

luminosity function from the latest APM survey of Loveday et al (1992). Thus dynamics of the 

Coma cluster implies the same value of the density parameter as we find from the independent 

analysis in this section. (Assuming that the mass-to-light ratio of galaxies varies according to the 

fundamental plane of ellipticals will result in < 10% correction, as will be discussed in the next 

section). 

To conclude the sections on the cluster correlation function and the pregalactic density field: 
1) We find that CDM models with any Q, h cannot fit simultaneously the data on the slope 
of the cluster correlation function, the dependence of its amplitude on richness and the APM 
angular correlation function. 2) We develop a method to determine the values of A(M)/As of the 
pregalactic density field directly from the data on the cluster correlation amplitude vs richness. 
The values determined in this way are independent of the numerical value of the biasing factor. 3) 
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The slope of the pregalactic density field, A(M)/A 8 vs the mass M, found in this way is consistent 
with that determined from APM. 4) Comparison of the amplitude of A(M)/A§ on a given linear 
scale r determined from the cluster correlation data with that determined from the APM survey 
fixes the value of $7. 5) The value of f2 determined in this way turns out ~0.25-0.3; the same 
value is implied by the dynamics of the Coma cluster. 6) If = 1 the amplitude of A(M)/A§ 
determined from the cluster correlation data does not pass through unity at 8/i _1 Mpc. 

7. Galaxy formation and high-z objects: constraining the small scale pregalactic 

density field 

In the previous sections we discussed constraints on the spectrum of the pregalactic density 
field from the present-day data. The smallest scales on which such constraints allow to probe the 
pregalactic density field come the cluster correlation amplitude - richness dependence and are 
around 5/i _1 Mpc (cf. Fig. 13). On smaller scales the spectrum of the pregalactic density field is 
constrained by observations of collapsed objects at high-z (Cavalieri and Szalay 1986; Efstathiou 
and Rees 1988 - hereafter ER88; Kashlinsky and Jones 1991; Kashlinsky 1993 - hereafter K93). 
Since the amplitude of the power spectrum is fixed by requirement that at z = it reproduce unity 
fluctuation in galaxy counts over a sphere of radius r$, for a given spectrum this normalization 
determines at what z the first objects of a given mass-scale would typically collapse. E.g. in 
the CDM models (eq.2) any increase in the amount of power on large scales would come at the 
expense of the power on small scales. Consequently, observations of collapsed objects at high 
redshifts when used in conjunction with the large-scale structure data can set strong constraints 
on CDM models (K93). In this section we discuss what the latest observational data imply. We 
also discuss the way the properties of the present day (elliptical) galaxies constrain the density 
fluctuations out of which they grew and collapsed. 

The data at high z which we use comes from observations of three types of objects at high 
redshifts: QSOs, galaxies and clusters of galaxies: 
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The latest grism surveys have been successful in finding quasars out to z ~ 5 (e.g. Schneider 
et al 1992,1994). There are now a couple dozen of quasars with z > 4 (e.g. Turner 1991 and 
references cited therein). This highest redshift comes from an optically selected quasar with 
z ~ 4.9 (Schneider et al 1991). There are suggestions that the true QSO abundance at high 
redshift may be significantly higher with most quasars hidden by dust obscuration (Fall and Pei 
1993). On the other hand, analysis of the grism surveys shows that the QSO comoving density 
peaks at z ~(2-3) and drops by a factor ~ 5 by redshift ~ 4.5 (Schmidt, Schneider and Gunn 
1991). This is also confirmed by data from the radio-loud quasars which are less likely to be 
affected by dust (Shaver et al 1996). Attempts have been made to use the QSO comoving density 
data at high z to test cosmological models (ER88; K93; Mahonen et al 1995). In order to do 
this one has to derive the total collapsed masses associated with high-redshift quasars and if, as 
is commonly assumed, QSOs hide an underlying galaxy, one has to further assume the efficiency 
with which galactic central nuclei lead to QSO phenomena (e.g. Begelman and Rees 1978, ER88). 
The lower bound on the mass comes from the Eddington luminosity limit on the total quasar 
power and the upper limit on the efficiency is < 1. If one assumes that a galaxy is hidden behind 
each quasar the total mass that collapsed at its redshift would be > 1O 12 M ; same number would 
be required by reasonable energy conversion efficiencies associated with the central engines (ER88; 
Turner 1991). In this case, low- CDM models would be the most difficult to reconcile with the 
QSO abundance at z ~(4.5-5), although Haenhelt and Rees (1993) argue that the required factors 
may be brought in agreement with CDM requirements at high z. While quasars can provide useful 
diagnostics of the various models, the uncertainties in their total mass, number densities as well 
as in the redshift of their collapse, rather than the redshift at which they are observed, make such 
constraints somewhat model dependent. 

Clusters of galaxies subtend linear scales comparable to r§ , where the rms density contrast 
today is unity. Therefore, existence of collapsed clusters of galaxies at high redshifts would imply 
a large initial density contrast on scale comparable to r$ and would be difficult to reproduce in 
models with £7 = 1 (K93). The data on the existence of clusters of galaxies at high z is only now 
becoming available, but it already indicates that galaxies may have assembled into clusters as rich 
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as Coma somewhere before z ~(l-2). Pascarelle et al (1996) reported a serendipitous discovery 
of a cluster or group of galaxies at z ~ 2.4 in the Hubble deep field. Francis et al (1996) report 
the discovery of a group of red, and therefore old, galaxies at roughly the same redshift. LeFevre 
et al (1996) discovered a cluster (or possibly group) of galaxies at z ~ 3.14. Dickinson (1993) 
reported observations indicating presence of rich galaxy clusters with a population of extremely 
red galaxies at z = 1.2. Recently, Jones et al (1997) detected a decrement in the temperature 
of the microwave background of 560/uK toward a pair of quasars at z ~ 3.8. Assuming this to 
be due to the Sunyaev-Zeldovich effect they searched for cluster of galaxies in that direction and 
concluded that the cluster must be at z > 1 with the total mass of ~ 1O 15 M (Saunders et al 
1997). Furthermore, such cluster mass and redshift are consistent with the quasar pair, observed 
to have very similar spectra and separated by ~ 10" , being gravitationally lensed images of the 
same quasar. Deltorn et al (1997) identified what is probably a massive cluster (~ 1O 15 M ) at 
z — 1 around a high redshift galaxy 3CR184 at the same redshift. The cluster is identified through 
both the gravitational arc near the central galaxy and by the excess of galaxies in the redshift 
distribution around 3CR184. All this suggests that clusters of galaxies may already have formed 
at z S> 1. Implications of the existence of such massive clusters at high redshifts for CDM models 
were discussed in K93. Fig.2c in K93 shows that such massive objects must be extremely rare 
fluctuations, >(7-10) standard deviations, in the density field of the CDM models. Most easily 
they could be explained by assuming low-0 and the pregalactic density field having power in 
excess of the CDM models at these scales. 

Galaxies at high redshifts are, probably, the most useful tool in constraining the power 
spectrum because in this case (at least) stellar mass can be fairly reliably estimated from colours 
and stellar populations. Until recently, only a handful of galaxies, all of them steep spectrum 
radio sources, have been found at high z, the most distant of which were 0902+34 and 4C41.17 at 
redshifts of 3.4 and 3.8 (Chambers et al 1990, Lilly 1988). Chambers and Chariot modelled their 
A-band photometry from Lilly (1988) with a rapid burst of star formation and concluded that 
the redshifts of formation must have been > 4. Eisenhardt and Dickinson (1992) have re-observed 
the 0902+34 galaxy at z = 3.8 in A-band and argue for its younger age. Eales et al (1993) have 
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discovered an optically luminous and massive radio galaxy 6C 1232+39 at z = 3.22. Lacy et al 
(1994) found a radio galaxy 8C1435+635 with a record z = 4.25. Djorgovski et al (1996) report 
the discovery of a galaxy responsible for a high-redshift damped Ly-a system at z = 3.15; from 
the velocity field they estimate the dynamical mass of the galaxy at > 2 x 10 11 h~ 1 M Q within 
15/t _1 Kpc. Hubble deep field observations of Steidel et al (1996) and Giavalisco et al (1996) find 
population of normal star forming galaxies at z > 3. Lowenthal et al (1997) have followed with 
spectroscopic observations at Keck of the Hubble deep field galaxies and with the 16 confirmed 
sources concluded that the comoving density of z > 3 galaxies is > 2 x 10~ 3 /i 3 Mpc -3 if Q = 0.1 
or > l(T 2 /i 3 Mpc- 2 if tt = 1. For comparison, the comoving density of L* galaxies today is 
~ 1.4 x 10~ 2 /i 3 Mpc (Loveday et al 1992). Trager et al (1997), using the Keck telescope, 
discovered four high-redshift galaxies: one at z = 3.35 and three at z ~ 4. Using the Ly-a break 
as redshift indicator Hu and McMahon (1996) reported the possible discovery of star-forming 
Ly-a emitting galaxies at z = 4.55. The above data indicate that galaxies must have been present 
already at z ~ 5. In an exciting recent discovery Franx et al (1997) identified a gravitationally 
lensed arc as galaxy at z = 4.92 and further find a companion to it with a radial velocity of only 
450 km/sec. Thus the data show that galaxies must already have been present in the Universe 
at redshifts beyond ~ 5. The existence of such early formation (collapse) of objects on scales 
~ 1O 12 M is indicative of small-scale power in the pregalactic density field in excess of that given 
by the low-fi CDM models required by the APM data (K93). 

Further indications of when galaxies must have collapsed come from recent observations 
of galaxies at smaller redshifts, but for which detailed spectroscopy and therefore accurate age 
determination are available. Hu and Ridgway (1994) found galaxy candidates whose colours 
were consistent with colours of ellipticals with old stellar populations shifted to z ~ 2.5 and 
have provided the first evidence that galaxies at high redshifts may already contain old stellar 
populations. A real breakthrough came with observations of the 53W091 galaxy at z = 1.55 
by Dunlop et al (1996, hereafter D96). D96 obtained high resolution spectroscopy of 53W091 
with the Keck telescope that revealed a very old stellar population in the galaxy. By fitting the 
spectrum with stellar models they determined its age to be £531^091 = 3.5Gyr with a very small 
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dispersion (0.5Gyr at 90% confidence level). This would rule out Einstein-de Sitter Universe for 
any H > 40km/sec/Mpc and for open or A-dominated flat models would place formation of 
53W091 at high redshifts z >(4-5). That this galaxy is unlikely to be either unique or unusual is 
evidenced by discovery of more such objects: e.g. 53W069 at z = 1.41 with age £53^069 = 4Gyr 
and 0.5Gyr uncertainty (Dey et al 1997). Francis, Woodgate and Danks (1997) report discovery 
two high redshift galaxies at z = 2.38 with old stellar population of > 0.5Gyr. 

Kashlinsky and Jimenez (1997, hereafter KJ97) discussed implications of the 53W091 data 
on low-f2 flat (f2 + A = 1 where A = CDM models. They pointed out that the redshift of 
formation of 53W091 required by its age would decrease with decreasing f2. But at the same, 
decreasing Q would in such models suppress the power on small scales thereby delaying galaxy 
collapse. KJ97 computed the total stellar mass of 53W091 by fitting the observed flux in all 
(V,J,H,K) bands with synthetic galaxy spectra based on the Miller-Scalo initial mass function for 
stars. The total mass in stars alone is given in Table 1 of KJ97 and is > 1O 12 M inside the aperture 
radius subtending (8-15)/i _1 Kpc at z = 1.55. The resultant stellar mass has little dependence 
on cosmological parameters or metallicity. KJ97 showed that in the flat A-dominated CDM 
cosmogonies 53W091 would have to be greater than ~ 5 standard deviations in the pregalactic 
density field. 

Further observational constraints on Q + A = 1 cosmologies come from the data on 
high-redshift type la supernovae. The latest measurements include 7 SNIa out to z — 0.46 and set 
limits of A < 0.35 at 68% or < 0.51 at 95% confidence levels (Perlmutter et al 1997). Similarly, the 
data on gravitational lenses implies A < 0.66 at 95% confidence level (Kochanek 1996). Together 
with the KJ97 results this makes flat A-dominated cosmogonies less attractive than open Universe 
cosmogonies. Therefore, in the remainder of this section we evaluate the limits on the pregalactic 
density field assuming A = 0, but most of the results can be relatively easily extended to A 7^ 0. 
For £1 + A = 1 Universe our main conclusions will not differ appreciably from the open Universe 
case. 

As discussed above the epoch of galaxy collapse constrains the pregalactic power spectrum 
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on galaxy scales and the data on the 53W091 galaxy places strong constraints on A-dominated 
flat Universe CDM models. In what follows we generalize the discussion of KJ97 to open 
Universe. Fig. 14 plots the redshift of formation of stellar population of 53W091 vs S7. Dotted 
lines correspond to the best estimate of its age from D96, t^woi = 3.5Gyr. Solid and dashed lines 
correspond to ±0.5Gyr uncertainty in £53^01 at 95% confidence level. Three lines of each type 
correspond to the Hubble constant values of h = 0.5,0.6 and 0.75. The age of 53W091 implies 
large values of the redshift of its formation, Zf or > 5, and agrees well with observations of other 
high redshift galaxies discussed earlier in this section. The numbers for another galaxy, 53W069 
(Dey et al 1997), would be similar, and both galaxies are inconsistent with the Q = 1 Universe for 
any reasonable value of the Hubble constant (h > 0.4). 

In order to estimate the likelihood of formation of objects like 53W091 in open CDM models 
we proceed in the manner outlined in K93. We start with density field in linear regime at some 
early epoch zf, in the formalism outlined in K93 and KJ97 the final results are independent of Z{. 
We denote with 5 co i the amplitude which density fluctuation had to have at Zi in order to collapse 
at redshift z. It can be evaluated from the spherical model approximation (e.g. Peebles 1980; 
Gott and Rees 1975; Narayan and White 1988; K93). The time evolution of a spherical shell that 
contained density contrast Si at Zi can be described by its "Friedman" equation: r 2 = A/r + C 
where A = H^Q, and — C = ii/g [|f2£j(l + Zi) + SI — 1] and r is expansion factor of the fluctuation. 
The turn-around time is given by t(r = 0) = ^A(-C) 3 ^ 2 and the collapse time is t co \ = 2t ta - The 
value of S co i is then determined by the condition that t co \ equals the cosmic time that elapsed 
between Zi and the redshift z when the fluctuation collapses: 



In order to normalize the density field at Z{ and also eliminate the dependence on Z{ we 
compute also the amplitude, As,*, which the fluctuation had to have at Zi in order to grow 
to 1/6 amplitude at z = 0. This can be done by noticing that conservation of mass requires 
1 + 5{z) = (1 + 5i)[(ri/r(z)] 3 . The fluctuation starts with Ag^ and must reach density contrast 




(34) 



of 1/b at z = 0. The time it takes is dt/f with r(z = 0),rj related by the above 
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mass-conservation expression. Hence the equation for A 8 j is given by: 

/■(l+l/6)-3 ^ r ^ /•*; dz 

For Zj>l and 5 co i, Ag^ <C 1 the ratio 5 co i/ Ag^ determined from eqs. (34), (35) is independent 
of Zi. We are interested in the pregalactic spectrum of density fluctuations and in linear regime 
fluctuations evolve independently of z. Hence, in what follows we omit the subscript i in the ratios 
(5/Ag as long as both amplitudes correspond to the pregalactic density field. 

Fig. 15 plots 5 co i/As vs 1 + z at which the fluctuation collapses for A = 0. Solid lines 
correspond to b = 1 and dotted lines to b = 2, but the solution of eq.(35) scales Ag,j oc 6 2 / 3 . The 
three lines of each type correspond to Cl = 0.1, 0.3 and 1 from bottom to top. One can see that the 
growth of density fluctuations slows down after 1 + z ~ Fig. 15 in conjunction with Zf or shown 
in Fig. 14 allows one to estimate the density field of 53W091 and the likelihood of its formation in 
a given cosmological model given the total mass of the galaxy. (For Q + A = 1 the ratio e) co //Ag is 
plotted in Fig.lc of K93; in that case the growth of fluctuations freezes at 1 + z ~ .) 

KJ97 estimated that the stellar mass of 53W091 must be at least 1O 12 M with little 
dependence on cosmology or star formation history. Thus in order to estimate the likelihood of 
53W091 in the open CDM models we adopt a (conservative) value for its total mass of 5 x 10 12 Mq. 
Fig. 16 illustrates the likelihood of 53W091 forming in the open CDM models. Thick lines in 
Fig. 16 show the value of the rms fluctuation, A(M), in units of Ag in the open CDM models on 
the (total) mass scale of 5 x 1O 12 M . At this mass range the CDM spectra have effective power 
slope index of ~ —2.5, so the ratio A(M)/Ag oc M~ 0A over the mass-scales relevant for 53W091. 
The three lines from top to bottom correspond to Hq =75,60 and 50 km/sec/Mpc. Thin lines 
show the bias-factor-independent values of S^wmi/ {^$b 2 ^) evaluated from eqs. (34), (35) at the 
redshift plotted in Fig. 14. Notation of the thin lines in Fig. 16 is the same as in Fig. 14. The values 
of the bias parameter for such models are determined by normalization to the COBE DMR data 
and are b > 1 (e.g. Stompor et al 1995). One can see from the figure that similarly to the flat, 
fl + A = 1, CDM models discussed in KJ97, the data on 53W091 would require it to be a very rare 
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(> 5 standard deviations) fluctuation in the density field of open CDM models. 

While the above arguments allow one to rule out or confirm certain cosmological models, they 
show only how rare the observed object is in a given model and thus do not allow to estimate the 
amount of power in the small scale pregalactic density field directly. Furthermore, estimating the 
number density of such objects requires assumptions about the gaussianness of the pregalactic 
density field. E.g. at a given number of standard deviations for gaussian models, such as CDM, 
collapsed objects would be less abundant than in models with non-gaussian density field such as 
strings (cf. Mahonen et al 1995). A different argument about galaxy formation epoch is required 
in order to estimate the small-scale pregalactic density field directly. 

We now move to estimating the epoch of galaxy formation from the present-day galaxies and 
then estimating the density field from which they formed directly. We assume that galactic haloes 
formed by dissipationless gravitational clustering (cf. White and Rees 1978). We assume further 
that the forming halo turns-around at radius r ta with little kinetic energy, i.e. its total energy 
per unit mass is E = —GM/rt a . The ensuing cold collapse will lead to violent relaxation during 
which the time-dependent gravitational potential redistributes energy and particle distribution 
(Lynden-Bell 1967), while the total mass of the halo remains conserved. Thus at the end of the 
violent relaxation the kinetic energy, T, and potential energy, V, will be related by the virial 
equilibrium relation: 2T + V = 0. I.e. the total energy per unit mass is E = T + V = —T = —\o 2 
where a is the one-dimensional halo velocity dispersion after virialisation. From the above 
the latter would be given by: a 2 = ^GM/r ta - We now assume that the halo turns-around at 
redshift z ta with an over-density 'jta(zta) with respect to the average density of the Universe, 

3H 2 

g^J7(l + zta) j at that epoch. Therefore for halo that has one-dimensional velocity dispersion a 
and, at least initially, had total mass M (which later could have been lost by e.g. tidal stripping) 
the redshift of its turn-around is given by: 

o2/3 1 

a 2 = —(GMH ) 2/3 in lt a(zta)}*(l + Zta) (36) 
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or in terms of numbers: 

nih^Mid + - 13.3 (^^) 2 (^^)" /3 (37) 

Once the redshift of the halo turn-around has been evaluated one can calculate the redshift of 
its collapse, z co i, by assuming the collapse to have taken twice the turn-around time (e.g. Gott 
and Rees 1975). As expected, at a given a the low mass haloes should have collapsed at higher 
redshifts. 

We turn to computing the turn-around overdensity Jtai 2 ) = Pta/p( z ta)- The turn-around 
time of a fluctuation is given by t ta = t(f = 0) = jjj^ [$ns-z-+n-i] 3 / 2 ' ^ n ^ e ^ e °^ ner hand, the 
linear scale at the turn-around is n a = rjO/(|$7c)jZj + 0—1), so the mean fluctuation density at 

3H 2 (-Cl5z+Q— l) 3 

the turn-around is given by pt a = g^c? ^ — t lT^ ■ Rewriting the right-hand-side of the above 

expression in terms of cosmic time at turn-around, tta, leads to: 

lta{z) = W) = (2) 0(1 TWW^W 2 m 

For = 1 one recovers from eq.(38) a well known result of Jta( z ) = (37r/4) 2 . For low values of O 
the density contrast at turn around is larger. Fig.17 plots jta vs z iov Q = 0.1 (dashed lines), 0.3 
(dotted lines) and 1. At 1 + z > Q^ 1 the values of jta converge to (37r/4) 2 . 

Eqs. (37), (38) allow one to estimate the redshift of galaxy collapse, given by t co i = 2tt a , 
for the halo of one-dimensional velocity dispersion a and mass M. Fig. 18 plots the value of z co \ 
vs O for galaxy with a = 200 km/sec. Solid line corresponds to M = W 12 h~ 1 M & , dotted to 
M = 2 x W 12 h- 1 M Q , and dashed line to M = 4 x lO 12 /i _1 M . If Q is low the redshifts at which 
such galaxy should have formed are in general agreement with the data discussed in this section. 
If, however, = 1 then massive galaxies would form at redshifts lower than the observations 
indicate. 

In fact, one can narrow down the range of z co i if one uses the observed properties of galaxies 
in order to compute the total halo mass and its velocity dispersion. In order to estimate the 
total halo mass at the time of its formation we proceed as follows: encouraged by the results that 
normalization to the observed properties of Coma gave in the previous section, we assume that 
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this cluster is representative of the global value of the mass-to-light ratios associated with galactic 
systems. In order to relate the initial mass of the halo to the luminosity of the galaxy it contains, 
we use the D n -a relation for elliptical galaxies (Dressier et al 1987; Djorgovski and Davis 1987). 
In obtaining our results we will thus be restricted to elliptical galaxies, but we will assume that 
the results are representative of the entire galaxy population. The fundamental plane relations of 
ellipticals are equivalent to the mass-to-light ratio of stellar populations of ellipticals scaling oc L B . 
From a detailed study of 37 ellipticals by van der Marel (1991) we adopt k = 0.35. Assuming 
further that the ratio of the dark-to-luminous matter is constant throughout the Universe and 
that luminous parts of galaxies formed by gas dissipation inside formed haloes (White and Rees 
1978) would allow to relate the total (initial) mass of formed haloes to the blue luminosity of the 
present day galaxies. (By "initial" we mean the mass at the time of the halo collapse; clustering 
processes would have later stripped parts, or all, of the halo material to form a continuous dark 
matter distribution of clusters and groups of galaxies.) Thus the mass of the halo at the time of 
its collapse is assumed to be given by: 

fM\ T(2-a) fL B \ 1+K fM\ f L B \ 1,35 , , 

M Halo = — \ „A -r) -H rr (39) 

V^B/Coma T{2 + K - a) \ ) \L B J Coma V ^* / 

The T-function factors in (39) come from requiring that (M/L) Co ma = J $(L)M(L)dL/ J (f>(L)LdL 
with $(L) = L- 1 $,(L/L,)- a exp(-L/L,) 

Observationally it is known that inner parts of elliptical galaxies are dominated by both 
stellar and halo potentials (e.g. Rix et al 1997) which would be consistent with ellipticals forming 
by dissipational collapse (Kashlinsky 1982). Thus central velocity dispersions of elliptical galaxies 
would not provide a good measurement of the halo potential well. The latter can be determined 
from dynamics in the outer parts of elliptical galaxies, where the only available measurements 
today come from X-ray data from hot gas coronae. The best measurements of temperature 
electron density profiles around ellipticals come from ASCA (Matsumoto et al 1997, Matsushita 
1997) and ROSAT (Davis and White 1996) satellites. Since the cooling time of the gas is of the 
order of Hq 1 the gas is likely to be hydrostatic. Then its temperature will be related to the halo 
velocity dispersion via the equation of hydrostatic equilibrium which for isothermal halo profile 
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. ^(^Boltzman rigasT) 



/inip 2 ^, where \x is molecular weight and m p is the proton mass (e.g. 



dr 



Awaki et al 1994). Temperature profiles (Matsushita 1997) indicate that the gas in outer parts 
can be assumed isothermal and that the gas has density profile n gas oc r~@ 9as with f3 gas ~(2.5-3). 
Then the halo velocity dispersion is related to the measured X-ray temperature via: 



Eqs. (37)-(40) along with the numbers in Fig. 15 allow now to determine the density fluctuation 
spectrum given photometric (Lb) and X-ray (T) data on galaxies. 

We have computed the values of 5(M)/A$ using eqs(37)-(40) for a sample of galaxies for 
which good X-ray and photometric data are available. The halo velocity dispersion was computed 
according to eq.(40) from the ROSAT (Davis and White 1996) and ASCA (Matsumoto et al 1997) 
observations. The data on the blue absolute luminosities were taken from Faber et al (1989) 
who determined the absolute luminosities using the fundamental plane relations. Only galaxies 
which appear in both the X-ray and Faber et al catalog were used. The ROSAT determined 
X-ray temperatures and gas profiles for galaxies have larger uncertainties than ASCA; hence of 
the ROSAT data listed in Table 1 of Davis and White (1996) we kept only those galaxies that 
have quoted uncertainties in T of less than 25%. The ASCA data available are listed in Table 3 
of Matsumoto et al (1997). Note that the atomic modelling of the spectra alone can introduce 
systematic uncertainties of ~ 20% (Matsushita 1997). The final dataset for which we found the 
accurate X-ray temperature measurements and photometry included some 20 galaxies. 

Fig. 19 plots the values of 6(M)/Ag vs the initial halo mass computed from eqs.(37)-(40) 
assuming a "universal" hot gas profile of f3 gas = 2.5 and fi = 0.6; the numbers in eq.(40) are not 
very sensitive to the value of (3 observed to lie between 2 and 3. The left box shows the results for 
£1 = 0.3, as required by the Coma cluster dynamics. The right box assumes the initial dark haloes 
to have the total mass-to-light ratio like Coma, but the total Q to be 1. Note that this prescription 
gives the actual amplitude of the density fluctuation, rather than its rms value A(M), and hence 
should have larger scatter. Additional scatter in Fig. 19 is introduced by the uncertainties in the 
current X-ray data and modelling and a further ~ 15% uncertainty in the individual distances 



i 




(40) 
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from using the fundamental plane relations. The slope in Fig. 19 is quite obvious and is in good 
agreement with the slope determined for large scales. The scatter is larger though, but is expected 
given the above observational uncertainties, the simplifying assumptions that were used in this 
section and the fact that the actual <5(M) rather than the rms, A(M) is shown. On the other 
hand, the X-ray temperatures do not correlate with the central velocity dispersion, and so the 
Faber-Jackson blue luminosity - central velocity dispersion correlation is unlikely to reproduce 
the trend in Fig.19. It is likely that the numbers in Fig. 19 reflect the slope and amplitude of the 
pregalactic density field on the corresponding mass scales. 



8. Reconstructed density field and cosmological paradigms 

The methods outlined in sections 3-7 allowed us to reconstruct pregalactic density field from 
various astronomical data on scales from ~ 1 to lOOft^Mpc. Fig. 20 summarizes the results in 
terms of the quantity: 

A(M) 
A 8 

with r being the linear scale containing mass M according to eq.(32). Left box corresponds to 
SI = 0.25, the value suggested by dynamics of the Coma cluster and galaxy luminosity function. 
Right box plots the reconstructed field for fi = 1, implied by inflationary prejudices. 

The open box in the figure shows the "normalization" point where the above quantity is unity 
by definition and irrespective of the bias factor. The three dotted lines show the BE93 fit to the 
APM data within one standard deviation uncertainty and the three lines correspond to the less 
accurate, but earlier, K92 fits to the data. The lines are drawn following the peculiar velocities 
analysis in sec. 4 suggesting the constancy of the bias parameter, b. They are plotted on scales > r$ 
where the APM data are less affected by non-linear gravitational effects. On scales > 100/i _1 Mpc 
the APM data probably are not useful for reliably probing the pregalactic density field. The 
rhombs show the values for A(M)/Ag reconstructed in Sees. 5, 6 from the data on the cluster 
correlation amplitude vs richness relation and normalizing the cluster richness-mass relation to 
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Coma. This reconstruction allows to determine the pregalactic density field on scales from ~ 5 to 
~ 20/i~ 1 Mpc. The two methods should give the same results where they reconstruct the density 
field on the same scales. For SI = 0.3 the density fields reconstructed from the two independent 
methods coincide. If O = 1 the discrepancy between the two amplitudes is significant, although 
the slopes coincide. Thus we conclude that O ~(0. 25-0.3) is in better agreement with the data and 
the pregalactic density field on scales 5/i _1 Mpc < r < 100/i -1 is likely to be close to that in the 
left box of Fig.20. 

The plus signs in Fig.20 show the density field reconstructed in Sec. 7 by using the fundamental 
plane relations and X-ray data for measuring halo velocity dispersion. The total halo masses at 
the time of their collapse/formation were normalized to reproduce the observed blue mass-to- light 
ratio of Coma. Plotted is the value of 8(M), not the rms A(M), in units of As; the numbers 
plotted depend on the amplitude of the bias parameter oc b 2 / 3 . Both the slope and amplitudes 
of the density field on galactic scales are consistent with extrapolation from the density field 
reconstructed on larger scales assuming a simple power law for P(k) with the slope of ~ —1. 
The slight excess in the amplitude (~ 30%) of the plus signs should not be regarded as worrying 
since what is plotted is the actual S(M), not its rms value. Furthermore, the galaxies used in the 
fundamental plane and X-ray studies are found predominantly in rich clusters and may indeed 
have formed out of higher peaks of the density field. Also the X-ray temperatures determined 
from ASCA and ROSAT observations for galactic X-ray emission may have systematic errors of 
~ 20% due to uncertainties in the atomic physics involved (Mushotzky, private communication). 
Such systematic errors may be sufficient to bring the amplitudes of the plus signs by the necessary 
~ 30%. In addition, many of the assumptions used here may be too simplifying: e.g assuming the 
same slope of the gas profile, P gas , or there may be deviations from the isothermal density profile, 
Phalo oc f~ 2 assumed here for all haloes. Given these uncertainties we find that the consistency in 
the density field plotted in Fig.20 from the various datasets and methods is quite reasonable. 

Fig. 21 juxtaposes the reconstructed density field with the predictions of the CDM models for 
the Harrison-Zeldovich (right box) and tilted spectra. Notation of the data points is the same 
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as in Fig.19. Only the data points for f2 = 0.3 (left box in Fig. 19) are shown because of the 
discrepancy between the cluster correlation data and the APM methods for Q = 1. CDM models 
are shown with dashed lines: thick dashed lines correspond to Vth = 0.5 and thin lines to Qh = 0.2 
for n = 1 and to 0.3 for n = 0.7. One can see that the models that fit the large-scale part of the 
pregalactic spectrum, miss the points on small scales by a significant factor. That is, for any value 
of the excess power, Qh, or n one would have to assume a significant scale dependence in the bias 
parameter b in order to fit the density field over the entire range of scales. 

The most economical conclusion to draw from the points in Fig.20 would be to assume that 
1) the Universe is open with O ~0.25-0.3, the value which is also suggested by the dynamics of the 
Coma cluster; and 2) the pregalactic power spectrum has constant power index of ~ — 1 all the way 
from ~ 30/i _1 Mpc to ~ l/i _1 Mpc. For primordial adiabatic fluctuations this could be achieved 
only if the smaller wavelengths which enter the horizon in the radiation dominated era have 
more power when compared with the scale- free primordial spectrum, Pi(k) oc k n . Computations 
show that this may be done with strings (Albrecht and Stebbins 1992) although the models are 
still not definite enough to predict a unique power spectrum. Another possibility could be the 
primeval-baryon-isocurvature (PBI) model where the shape of the spectrum of the density field 
could be further evolving after recombination due to reionisation and other effects (Peebles 1987). 
Or one can assume that the primordial power spectrum in the CDM models was not scale free 
and contained extra power on small scales. However, the latter may lack theoretical motivation if, 
as we have argued on the basis of cluster correlation and galaxy data the Universe is open, and 
CDM model predictions are tied to and/or motivated by inflationary scenario. 

Fig.20 plots the pregalactic density field in units of A§ and thus independently of the latter's 
value. In order to assign a numerical value to A(M) at some early epoch, say the epoch of 
recombination, we have to compute the value of As at that time. Assuming purely gravitational 
evolution, i.e. ignoring pressure forces, Thompson drag and other effects which can become 
important in some models, would give at z = 1000 the values of Ag = (0.55, 1.2,2.8) x 10~ 3 6 2 / 3 
for = (1,0.3,0.1) respectively. The value of b can be determined from normalization to the 



-50- 



microwave background anisotropy measurements. We omit this since it would require assumptions 
about the nature of the density field, dark matter, reionisation history, etc. 

9. Conclusions 

The results of the paper can be summarized as follows: 

1) Our analysis of the APM survey results on w{6) in six 0.5 magnitude wide slices shows 
that the various evolutionary corrections are not significant out to the limits of the catalog 
at bj ~ 20.5. This allows to make accurate fits of the various spectra to the APM data. Of 
these, the BE93 fit for the power spectrum of the density field traced by galaxies gives the best 
approximation to all six slices. If light traces mass, this fixes the pregalactic density field on scales 
10/t _1 Mpc < r < 100/i -1 Mpc. We then choose the value of the angular scale at which the angular 
correlation function is ~ 3.3 times the systematic error from the plate gradients in all six slices to 
test the requirements the data sets on CDM models. We find that in order to fit the APM data, 
CDM models would require Uh ~ 0.2 if n = 1 or Qh ~ 0.3 if n = 0.7. 

2) The velocity field implied by the galaxy correlation function determined from the APM 
data is consistent with the Great Attractor peculiar velocities and the bias parameter which is 
constant over the range of scales probed by both datasets, 10/t _1 Mpc < r < 60/t _1 Mpc. This 
suggests that pregalactic density field is proportional to the power spectrum of galaxy clustering 
determined from the APM on scales > r§. Furthermore, if Q 6 /b = 1 the amplitude of the dot 
velocity correlation function at zero lag would be significantly larger, ~ (1300km/sec) 2 , than the 
values indicated locally and by the central velocity dispersions of typical galaxy systems. This 
suggests that the velocity data are in agreement with open Universe. Furthermore, we show that 
the quantity constructed from the the velocity correlation tensor, ^/S(r) — n(r), can be used to 
determine fi a6 /6. Its amplitude, predicted by both the APM and and the (local) CfA data on 
£(r), is ~ 400^^- km/sec on scales < 40/?T 1 Mpc and should be detectable in the current velocity 
surveys if Q°- 6 /b = 1. 
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3) We calculate in detail the amplification in the cluster correlation function for clusters of 
various masses due to gravitational clustering assuming only that clusters of galaxies should be 
identified with regions with turn-around time less than the age of the Universe. We use the data 
on both the slope of the cluster correlation function for Abell clusters of richness 1Z > 1, and on 
the dependence of the cluster correlation amplitude on the cluster richness/mass. Cluster masses 
for given 1Z are normalized to Coma. We find that, for any value of Qh, CDM models cannot 
simultaneously fit data on the cluster correlation slope and amplitude, its amplitude-richness 
dependence and the APM catalog data on w(9). We show that by using the data on the underlying 
galaxy correlation function, £(r) at some fixed r (we chose r = 25/i~ 1 Mpc), enables one to invert 
the amplitude-richness data to obtain the values of the rms fluctuations, A(M)/Ag, for the 
pregalactic density field on cluster mass-scales. Applying the method to the data we obtain the 
pregalactic density field whose slope is in excellent agreement with the APM field on the same 
scales. Requiring then that the two amplitudes coincide fixes f2 to be 0.3 independently of the 
bias factor; the same value that is implied by the dynamics of the Coma cluster. 

4) We then use the data on objects at high redshift to constrain the small-scale part of the 
pregalactic density field. We argue that data on the existence of clusters of galaxies at redshifts 
>(l-2), observations of galaxies at z ~ 5 and the recent determination of galaxy ages from Keck 
observations for a number of high-z galaxies imply more power on small scales than predicted 
in CDM models normalized to the large-scale galaxy correlation data. We then reconstruct the 
pregalactic density field out of which modern-day (elliptical) galaxies have formed. To do this we 
use the data on blue absolute luminosities of Faber et al (1989) and assume the fundamental plane 
relations to determine the initial halo mass (normalized further to reproduce the mass-to-light 
ratio of the Coma cluster) and the X-ray data from ROSAT and ASCA satellites to determine the 
halo velocity dispersion. Despite the simplifying assumptions we recover a density field which is in 
good agreement with the simple extrapolation of the reconstructed density field from larger scales. 

5) At the end of these steps we reconstruct the pregalactic density field over linear scales 
encompassing two orders of magnitude, l/i _1 Mpc < r < 100/i _1 Mpc. The recovered field is 
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difficult to fit with CDM models and constant bias parameter over this range of scales. We argue 
that the most economical explanation of the density field on all scales would be to assume that 
the results in Fig. 20 represent the pregalactic density field with constant bias parameter b ~ 1 
and 12 ~ 0.3 as required by dynamics of galaxy systems and the consistency between the X-ray 
observations of rich clusters and Big-Bang nucleosynthesis (White et al 1993) 
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Figure captions 

Fig.l: Selection function for the six 0.5 magnitude wide slices of the APM data. The lines 
are computed using the if -correction as Sm = Kz with K = 3, but the dependence on K is weak 
for K < 4. Solid lines correspond to Q = 0.1 and dotted to Q = 1. The six curves of each type 
correspond to the magnitude range of each of the APM slices shown on top of the boxes in Fig.2. 

Fig. 2: Open triangles show the APM data in the narrow magnitude slices with the magnitude 
ranges shown on top of the boxes. Thick dashed lines are for Q = 1 and a power law £(r), which 
should give a good fit to the data at small angular scales. Thick dotted lines, that mostly merge 
with the thick dashed lines, is the power law fit for Q = 0.1. The upper thick lines of each type 
are for the clustering pattern stable in physical coordinates, i.e. £(r;z) = (r /r*)~ L7 (1 + z)~ iri . 
Lower thick lines correspond to clustering pattern which is stable in comoving coordinates: 
£(r; z) = (r /r*)~ L7 '(1 + z)~ 3 . Three solid lines correspond to the best BE93 fit to the APM power 
spectrum and ± one standard deviation uncertainty. For the BE93 lines the power spectrum was 
taken to evolve oc (1 + z) -3 . The numbers are shown K = 3. 

Fig.3: The angular scale, 9(w), at which the projected angular correlation reaches a given 
value, w, is plotted vs the mean blue magnitude of the APM six slices. Left boxes correspond 
to w = 2 x 10 -3 which is just above the systematic error, w pp ~ 1.5 x 10 -3 . The right boxes 
are for w = 5 x 10~ 3 . Thick plus signs correspond to the APM data plotted in Fig.2. The data 
is compared with CDM predictions: upper boxes show the values of 9(w) for n = 1 and lower 
ones for the tilted CDM model with n = 0.7. x signs correspond to Qh = 0.5, open squares to 
Qh = 0.3, triangles to tth = 0.2, rhombs to Qh = 0.1 and asterisks to Qh = 0.05. Only CDM 
models which match the observed values would be consistent with APM data. 

Fig.4: The values of the rms density contrast, A(M), in units of Ag are plotted for the various 
fits to the APM data. Dotted lines show the BE93 fit with the one standard deviation uncertainty. 
Solid lines correspond to the K92 fit with ko = 35, 40, 50/i _1 Mpc from bottom up. Thick dashed 
line corresponds to CDM models with Qh = 0.2, n = 1 and Qh = 0.3, n = 0.7. 
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Fig.5: The values of the "typical velocity" profile derived from eq.(ll) using the data on the 
galaxy correlation function, £(r), is plotted for Q os /b = 1. Dashed line corresponds to the velocity 
field from the CfA data on £(r). Solid lines correspond to £(r) given by the BE93 fit to the APM 
data with one standard deviation uncertainty. 

Fig.6: The values of the transverse minus parallel velocity correlations according to the data 
on £(r). Notation is the same as in Fig. 5. 

Fig.7: Cluster-cluster correlation function for CDM models with n = 1 and various Q,h on 
scales 20/i _1 Mpc < r < 100/i -1 Mpc. Thin solid lines show the underlying correlation function 
in the CDM model specified on top of each box. Dotted lines shows the amplified correlation 
function for clusters of richness N = 5, dashes to N = 10, dashed-dotted lines to N = 50, 
and the dashed-double-dotted lines to N = 100. For comparison, the richness class 1Z > 1 
clusters, such as Coma, have N = 100. Thick solid line corresponds to the approximation 
£ cc = 300(r/l/i- 1 Mpc)- L8 which Bahcall and Soneira (1983) find for K > 1 clusters. 

Fig.8: Same as Fig.7 only for tilted CDM models with n = 0.7. Same notation as in Fig.7. 

Fig. 9: Amplification of the correlation function for clusters A = £mm(t")/£(t) evaluated for 
CDM models at r = 25/i~ 1 Mpc is plotted vs cluster richness. Triangles correspond to the data 
points from Fig.3 of Bahcall and West (1992). Solid lines are for CDM models with VLh = 0.1, 
dotted for Clh = 0.2, dashed for Qh = 0.3, and the dashed-dotted lines correspond to flh = 0.5. 
Thin lines of each type are for h = 1 and thick lines are for h = 0.5. 

Fig. 10: Correlation function £(r) according to the various fits to APM is plotted vs r. 
Dotted line correspond to CDM model with n = 1 and Qh = 0.2 or n = 0.7 and Qh = 0.3 
(the two practically coincide). Dotted lines correspond to the BE93 fit with one standard 
deviation uncertainty. Three solid lines correspond to the K92 fit with the transition scale 
k~ 1 = 35,40,50/i- 1 Mpc. 

Fig. 11: The amplification factor A(x) is plotted according to eq.(29). Solid line corresponds 
to the underlying correlation £ = 0.07 and the two dotted lines show the results of, say, 25% 
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uncertainty in £. 

Fig.12: The spectrum of the pregalactic density field, A(M)/Ag, obtained from the data 
on the cluster correlation amplitude - richness dependence. The upper horizontal axis shows the 
richness corresponding to the particular value of A(M)/A 8 . The lower horizontal axis shows the 
values of M computed from normalization to Coma (eq.31). 

Fig. 13: Plots A(M)/A 8 from Fig.12 vs linear scale (eq.32). Left box: thick plus signs 
correspond to 17 = 1, thin plus signs to 17 = 0.1. Open box corresponds to A(M)/A§ = 1 at 
r = 8/i _1 Mpc; the numbers miss this point if 17 = 1. Right box: dotted lines correspond to the 
BE93 fit to the APM data with one standard deviation uncertainty. Solid lines are for K92 fit 
with transition scale Uq 1 = 35, 40, 50/i~ 1 Mpc. Rhombs correspond to the data from Fig.12 with 
17 ~ 0.25 when it would coincide with the independently deduced lines from the APM data. Other 
signs are in the same notation as in Fig.12. 

Fig.14: The redshift of formation of 53W091 according to the age estimate t age = 3.5 ±0.5 Gyr 
of D96 is plotted vs 17 for A = 0. Dotted lines correspond to t age = 3.5Gyr, solid to t age = 3Gyr 
and dashes to t age = 4Gyr. Three lines of each type corresponds to h = 0.5, 0.6, 0.75 from bottom 
up. 

Fig. 15: Values of 5 co i/Ag at early epoch, Zi, vs redshift of collapse z are plotted for A = 0. 
Solid lines correspond to A§ evolving to 1 today, or b = 1; dotted lines are for 6 = 2. Three lines 
of each type correspond to 17 = 1, 0.3, 0.1 from up down. 

Fig. 16: The values of (5 co ;/Ag implied by the redshift of collapse of 53W091 plotted in 
Fig.14 are shown vs 17 for zero cosmological constant Universe. Three thick solid lines show 
the predictions of this ratio in CDM models with the value of 17 shown on the horizontal axis. 
The total mass of 53W091 there was taken to be 5 x 1O 12 M (cf.KJ97). The three thick lines 
correspond to h = 0.5,0.6,0.75 from bottom up. 

Fig. 17: The overdensity of the fluctuation at turn-around plotted vs the turn- around redshift. 
Solid line is for 17 = 1, dotted for 17 = 0.33 and the dashed line is for 17 = 0.1. 
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Fig. 18: The redshift of collapse (formation) of galaxy with the halo dispersion a = 200km/sec 
plotted vs Q according to eq.(37). The total (initial halo) mass is taken to be 1O 12 /i _1 M (solid 
line), 2 x HP/i^Mq (dotted line) and 4 x 10 12 fe -1 M Q (dashed line). 

Fig.19: Density fluctuations <5(M)/A 8 computed according to eqs.(37)-(40) vs the mass scale 
they contain. The numbers were computed normalizing the total mass-to-light ratio of galaxies to 
Coma at M/L = 362/i. Left box assumes fl = 0.3, corresponding to the value implied by Coma 
dynamics, and the right box corresponds to Q = 1. 

Fig.20: Pregalactic density field, A(M)/Ag reconstructed over the entire range of 
l/i _1 Mpc < r < 100/i _1 Mpc. Left box corresponds to Q = 0.3, right box to = 1. Open square 
corresponds to the value of 1 at r§. The lines from the APM data are plotted only for linear scales 
r > r§. Solid lines are the K92 fits to the APM data with the transition scale 1 = 35,40, 50h^ 1 
from bottom up. The three dotted lines correspond to the BE93 fit with one standard deviation 
uncertainty. Rhombs show the values obtained from inverting the data on cluster correlation 
amplitude vs richness according to eq.(29). The rhombs coincide with the APM data only if 
SI ~ 0.3. Plus signs show the values reconstructed according to eqs.(37)-(40) from the data on 
the fundamental plane of galaxies and the ROSAT and ASCA measurements of the halo velocity 
dispersion. 

Fig.21: The results from the left box of Fig.20 are juxtaposed with CDM models. In the left 
box the dashed lines show the CDM density field with n = 1 and Qh = 0.5 (thick dashed line) and 
Qh = 0.2. The right box shows the tilted CDM models with n = 0.7 and flh = 0.5 (thick dashed 
line) and flh = 0.3. Other symbols have the same notation as in Fig.20. 
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Table 1. Data from Bertschinger et al (1990) 

r = r = 40/i _1 Mpc r = 60/i -1 Mpc 

T/(r) (km/sec) 457 ± 61 388 ± 67 327 ± 82 

(L,B) (156 ±7; -19 ±12) (177 ± 9; -15 ± 17) (194 ± 13; 5 ± 26) 
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Table 2. y/v{0) in km/sec computed according to eq.(12) for the various fits to £(r). 
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